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NOTATION 


a  meridional  radius  of  curvature  for  a  toroidal  shell  with  circular 

cross  section  (see  Fig.  VI  -  1) 

b  distance  between  the  center  of  the  cross  section  and  the  axis  of  a 

toroidal  shell  (see  Fig.  VI  -  1) 

C  =  cos  nub 

m  T 

D  =  — ,  extensional  stiffness 

1  -  v 

E  Young's  modulus 

h  thickness  of  shell 

k  =  1/12  (h/a)2 

Eh3 

K  =  - »  bending  stiffness 

12(1  -  v) 


m  Fourier  index 

M  number  of  terms  used  in  series  expansions  for  the  displacement 

components 

n  number  of  circumferential  waves  in  buckle  pattern 


N0  •  N0  '  N0£ 


N<*>0  1  Q<*>  •  Q0 

, 

M<pe  •  Mo<t> 
N<p0  ’  Neo  ’ 
N^>0O’  N0^O 


p 


>  incremental  stress  resultants  (see  Fig.  Ill  -  2) 


prebuckling  stress  resultants  (see  Fig.  IV  -  1) 
external  pressure  loading  for  a  toroidal  shell 
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p  pn  ,  p  components  of  the  applied  loading  per  unit  area  of  the  shell's 
o  u  z 

midale  surface  (see  Fig.  IV  -  1) 
r  radius  of  a  parallel  circle  (see  Fig.  Ill  -  1) 

r,  ,  r9  meridional  and  circumferential  radii  of  curvature  (see  Fig.  Ill  --  1) 

S  =  sin  mii 

m 

t  time 

u  (</>),  v  (9)  , 

n  n  displacement  functions 

wn(4>) 

u  ,  v  ,  w  circumferential,  meridional,  and  radial  displacement  components 

of  a  point  in  the  shell's  middle  surface  (see  Fig.  Ill  -  3) 

U  ,  V  ,  W  Fourier  coefficients  [see  Eci.  (VI -35)1 

U  ,  V  ,  W  Fourier  coefficients  [see  Eqs.  (VI  —  20)  J 

m  m  m  1  M  ’ 

Xj  =  sin  9 

yx  ,  y2  -  cos  0  ,  cos  29 

z  radial  coordinate  (see  Fig.  Ill  -  2) 

a  -  b/a 


€<p’  €0'  y9 e 

€<f>'  €0'  ~ytp 0 

K0 ’ K0’ K9 0 
A 


e ,  9 


Kronecker  delta  [see  Eq.  (VI-25) J 
defined  in  Eq.  (IV  -  4) 
defined  in  Eq.  (VI  -  27) 

strains  at  a  distance  z  from  the  middle  surface 
middle  surface  strains 
curvature  changes 

-  pa/Eh  ,  nondimensional  pressure  parameter  for  a  toroidal  shell 
circumferential  and  meridional  coordinates  of  a  point  on  the 
shell’s  middle  surface  (see  Fig.  Ill  -  1) 
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Poisson’s  ratio 

mass  per  unit  area  of  the  shell’s  middle  surface 

-  p  -  ^  ,  meridional  coordinate  for  a  toroidal  shell  (see  Fig.  VI  -  1) 
.  1.  2 

=  ,  nondimensional  frequency  parameter  for  a  toroidal  shell 

-r  for  stability  of  a  toroidal  shell 
<  or 

circular  frequency  for  vibration  of  a  toroidal  shell 

» 

angles  of  rotation  of  the  normal  to  the  shell's  middle  surface 
(see  Fig.  Ill  -  3) 

angles  of  rotation  around  the  normal  to  the  shell's  middle  surface 
(see  Fig.  Ill  -  3  and  Eqs.  (III-10)  and  (III-l  1)  ] 
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INTRODUCTION 


The  first  solution  in  the  field  of  buckling  of  thin  shells  was  given  in  Lorenz's 
paper  on  axially  symmetric  buckling  of  axially  compressed  cylinders  (Ref.  1). 

In  1932,  Flugge  (Ref.  2)  developed  a  general  theory  for  buckling  of  cylinders  and 
presented  numerical  results  for  simply  supported  cylinders  under  various  loading 

| 

conditions.  For  the  case  of  the  uniformly  loaded  cylinder,  Flugge  made  the 

t; 

V 

usual  assumption  that  the  prebuckling  stresses  could  be  approximated  by  a 

1: 

homogeneous  membrane  state  of  stress.  Therefore,  the  stability  equations  con- 

» 

tain  constant  coefficients,  and  exact  solutions  can  readily  be  obtained  for  isotropic 

{ 

or  orthotropic  cylinders  with  arbitrary  boundary  conditions .  However,  when  the 

i 

applied  loading  is  nonuniform  the  differential  equations  governing  the  stability  of 

i 

the  cylinder  will  have  variable  coefficients.  Consequently  it  becomes  considerably 
more  difficult  to  obtain  solutions;  for  example,  see  Flugge' s  analysis  of  a 
nonuniformly  compressed  cylinder  (Ref.  2). 

For  shells  other  than  cylindrical,  the  stability  equations  contain  variable 
coefficients.  Most  of  the  work  done  on  noncylindrical  shells  has  been  devoted 
to  spheres  (Refs.  3  through  5)  and  cones  (Refs.  6  and  7).  Relatively  little  atten¬ 
tion  has  been  devoted  to  shells  with  variable  Gaussian  curvature.  Mushtari  and 
Galimov  (Ref.  8),  using  shallow  shell  equations,  presented  a  simple  formula  for 
the  critical  normal  pressure  of  an  ellipsoidal  shell.  Their  analysis,  however, 
appears  to  be  greatly  oversimplified.  Machnig  (Refs.  9  and  10)  investigated  the 
stability  of  a  torus  subject  to  uniform  external  pressure.  In  his  first  paper  (Ref.  9). 

Machnig  studied  both  axially  symmetric  and  asymmetric  buckling  modes  and 
concluded  that  ihe  former  gives  the  smallest  critical  pressure.  A  perturbation 
technique  was  used  to  solve  a  system  of  partial  differential  equations  governing 
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the  asymmetric  buckling  mode.  Of  course,  it  should  be  possible  to  separate  the 
space  variables  in  the  stability  equations  for  a  complete  shell  of  revolution  sub¬ 
ject  to  axially  symmetric  loads  and  thereby  obtain  a  system  of  ordinary  differ¬ 
ential  equations .  Apparently,  Machnig  had  to  contend  with  partial  differential 
equations  instead  of  ordinary  differential  equations  since  it  appears  that  it  is  not 
possible  to  separate  the  space  variables  in  his  equations  for  the  asymmetric 
buckling  mode.  In  his  more  recent  paper  (Ref.  10),  Machnig  considers  only  the 
axially  symmetric  mode.  A  review  of  that  paper  was  given  by  Koiter  (Ref.  11). 
The  reviewer  indicates  that  asymmetric  buckling  modes  may  well  result  in 
smaller  buckling  pressures  for  some  values  of  the  shell's  geometric  parameters 
and  that  the  power  series  expansion  must  break  down  for  toroidal  shells  with 
small  values  of  b/a  (see  Fig.  VI-1*  for  notation).  In  spite  of  these  critical 
remarks,  Koiter  concludes  that  Machnig' s  paper  must  be  regarded  as  a 
significant  first  step  in  the  solution  for  the  buckling  of  a  torus. 

In  this  work,  the  stability  of  a  general  shell  of  revolution  subject  to  arbitrary 
loads  will  be  investigated.  First,  the  elastic  law  which  relates  the  incremental 
stress  resultants  to  the  incremental  displacement  components  will  be  derived. 
This  law,  derived  from  elementary  considerations,  turns  out  to  be  the  same  as 
the  elastic  law  derived  by  Reissner  (Ref.  12)  who  used  the  methods  of  differential 
geometry.  Next,  the  equations  of  equilibrium  will  be  applied  to  a  differential 
element  of  the  deformed  shell.  The  resulting  partial  differential  equations  are 
linear  and  homogeneous  in  the  incremental  quantities,  and  the  specification  of 
linear  and  homogeneous  boundary  conditions  results  in  an  eigenvalue  problem. 


♦Combined  Roman-Arabic  numbers  designate  cross-chapter  references  of 
figures  or  equations. 
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Thereafter,  only  complete  shells  of  revolution  under  axially  symmetric  loads 
will  be  considered.  This  means  that  the  coefficients  in  all  equations  are  inde¬ 
pendent  of  the  circumferential  coordinate  9  and  that  all  incremental  quantities 
are  periodic  functions  of  6  (see  Fig.  HI-1  for  notation).  Therefore,  it  is  pos¬ 
sible  to  express  all  incremental  quantities  as  Fourier  series  in  9  and  replace 
the  partial  differential  equations  by  ordinary  differential  equations.  Then,  the 
stress  resultants  will  be  eliminated  from  the  equilibrium  equations,  and  three 
ordinary  differential  equations  for  the  three  incremental  displacement  components 
will  be  obtained. 

Next,  by  specialization  of  these  equations,  the  ones  governing  the  stability 
of  a  toroidal  shell  subject  to  a  uniform  external  pressure  are  obtained.  It  turns 
out,  as  might  be  expected,  that  two  types  of  asymmetric  buckling  modes  exist: 
one  which  is  symmetric  with  respect  to  the  plane  4  -  0,  n  (see  Fig.  Vl-1)  and 
one  which  is  antimetric  with  respect  to  this  plane.  This  will  be  evident  when 
series  expansions  are  used  for  the  incremental  displacement  components ;  as  a 
result,  two  uncoupled  systems  of  linear  homogeneous  algebraic  equations  are 
obtained  for  the  free  constants  in  the  series  expansions .  A  matrix  iteration 
technique  is  employed  to  obtain  the  lowest  eigenvalue  and  the  corresponding 
eigenvector  for  each  system.  After  the  eigenvector  has  been  computed,  the  mode 
shapes  are  given.  The  matrix  iteration  method  is  programmed  for  the  IBM  7094 
digital  computer  and  numerical  results  are  presented  in  the  form  of  design 
curves  which  give  nondimensional  buckling  pressures  for  a  wide  range  of  the 
shell's  geometric  parameters.  These  results  are  compared  with  available  test 
results  as  well  as  with  Machnig's  results. 

Finally,  equations  for  the  free  vibrations  of  a  pre-stressed  shell  of  revolu¬ 
tion  will  be  derived  and  then  specialized  for  the  toroidal  shell  subject  to  external 
or  internal  pressure. 
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II 


FUNDAMENTAL  CONCEPTS  AND  ASSUMPTIONS 


1.  Concept  of  Buckling 

In  order  to  establish  basic  notions  and  at  the  same  time  to  introduce  a 
certain  terminology,  we  will  disc.ss  the  concept  of  bifurcation  buckling  in  this 
section.  The  ideas  expressed  here  may  be  found  in  many  works;  e.g.,  Refs. 

13  through  18. 

Let  us  consider  a  conservative  mechanical  system  subjected  to  applied 
loads.  The  distribution  and  direction  of  the  loads  are  assumed  to  be  known. 

The  magnitude  of  the  loads  is  assumed  to  be  characterized  by  a  single  non¬ 
negative  scalar  X ,  which  we  call  the  load  parameter .  When  X  is  zero,  all 
the  applied  loads  vanish  and  the  system  is  undeformed.  For  other  valuer  of  X, 
the  system  will  deform  and  develop  stresses  in  order  to  be  in  equilibrium  under 
the  applied  loads.  Let  the  load  parameter  increase  monotonically  from  its 
initial  zero  value.  Let  us  assume  that  for  sufficiently  small  values  of  X,  there 
exists  a  unique  solution  to  the  equilibrium  problem  for  the  system.  The  equi¬ 
librium  configuration  given  by  this  solution  is  called  the  basic  state  or  prebuckled 
state .  It  may  happen  that  for  some  value  of  X ,  say  X  =  X* ,  there  exists  another 
solution  or  equilibrium  configuration,  called  the  buckled  state,  which  is  infinitesi¬ 
mally  close  to  the  basic  stale.  If  such  a  value  X*  exists,  then  it  is  called  the 
critical  value  of  the  load  parameter.  The  objective  of  the  buckling  analysis  is 
to  determine  X*  .  The  existence  of  two  adjacent  equilibrium  configurations  for 
the  same  value  of  the  load  parameter  means  that  there  is  a  bifurcation  of  the 
basic  state.  The  stresses  and  displacements  of  the  prebuckled  state  are  called 
the  prebuckling  stresses  and  prebuckling  displacements.  At  the  critical  value 
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of  the  load  parameter,  the  differences  between  the  stresses  and  displacements 
in  the  prebuckled  and  buckled  states  are  called  the  additional  or  incremental 
stresses  and  displacements. 

In  the  present  analysis,  the  mechanical  system  is  a  shell  of  revolution;  the 
applied  loads  are  the  conservative  external  pressures  p0(X) ,  p^(X) ,  PZ(X) ; 
the  basic  state  is  assumed  to  be  a  membrane  state  of  stress ;  the  prebackling 
stress  resultants  are  denoted  by  N0Q,  N^Q,  ;  and  the  incre¬ 
mental  quantities  are  denoted  by  u,  v,  w,  N0,  ,  N^,  Q0,  , 


M0’  M0’  M<£0* 

In  the  mathematical  formulation  of  the  stability  problem,  we  assume  that 
the  shell  is  in  a  buckled  state  and  we  seek  the  smallest  value  of  the  load  param¬ 
eter  for  which  all  the  basic  equations  can  be  satisfied.  The  basic  equations  in 
shell  theory,  or  in  any  branch  of  solid  mechanics,  are  the  equations  of  equilib¬ 
rium,  the  equations  for  kinematics,  and  the  constitutive  equations  (Hooke's  law 
for  a  linearly  elastic  material).  Before  presenting  these  equations,  we  list  in 
the  next  section  the  assumptions  and  limitations  used  in  the  analysis. 

2.  Assumptions  and  Limitations 

(a)  The  shell  is  made  from  an  isotropic  and  homogeneous  material  which 
obeys  Hooke's  law  (linearly  elastic  material). 

(b)  The  thickness  of  the  shell  is  constant. 

(c)  The  thickness  of  the  shell  is  small  in  comparison  with  the  radii  of 
curvature  of  the  middle  surface  (thin  shell  approximation). 

(d)  Kirchhoff-Love  hypothesis 

(i)  A  straight  line  normal  to  the  middle  surface  before  deformation 
remains  straight  and  normal  to  the  middle  surface  after 
deformation  and  retains  its  original  length. 
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(ii)  The  normal  stresses  acting  on  surfaces  parallel  to  the  middle 
surface  are  small  in  comparison  with  other  stresses  and  may 
be  neglected  in  the  stress -strain  relations. 

(e)  All  incremental  quantities  are  infinitesimal. 

(f)  The  incremental  strains  are  small  in  comparison  with  the  incremental 
rotations . 

(g)  The  basic  state  may  be  approximated  by  a  membrane  state  of  stress. 
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Ill 


THE  ELASTIC  LAW 

The  relations  between  the  incremental  stress  resultants  and  the  incremental 
displacement  components  of  the  middle  surface  will  be  presented  in  this  chapter. 
These  relations,  which  represent  the  elastic  law  for  the  shell,  are  obtained 
through  combination  of  the  equations  of  kinematics  and  the  constitutive  equations. 
The  elastic  law  for  the  shell  buckling  problem  considered  here  must  be  identical 
to  that  for  the  problem  of  infinitesimal  bending  of  a  shell  of  revolution  subject  to 
external  loads.  Therefore,  we  could  use  the  elastic  law  derived  by  Flugge 
(Ref.  13)  in  his  analysis  of  the  linear  bending  of  a  shell  of  revolution  under 
external  loads.  However,  instead  of  proceeding  directly  to  Flugge' s  final  results, 
we  will  closely  follow  his  derivation  of  the  elastic  law  in  order  to  see  what  assump¬ 
tions  are  made  and  at  what  stage  they  are  made.  The  elastic  law  derived  here  will 

•• 

differ  from  the  one  derived  by  Flugge  because  we  will  deviate  somewhat  from 
his  derivation. 

1 .  Geometry  of  a  Shell  of  Revolution 

The  middle  surface  of  a  shell  of  revolution  is  shown  in  Fig.  la.  The  equa¬ 
tion  defining  the  meridian,  y  =  y(r)  or  y  =  y (d>),  is  assumed  to  be  given.  The 
lines  of  principal  curvature  are  the  meridians  (0  =  const.)  and  parallel  circles 
(rt  =  const.) .  A  point  P  on  the  middle  surface  of  the  shell  is  determined  by 
the  intersection  of  a  meridian  and  a  parallel  circle.  Accordingly,  we  rse  the 
angles  0  and  0  as  the  curvilinear  coordinates  of  the  middle  surface  of  the 
shell. 
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The  principal  radii  of  curvature  r  ^  and  r2  are  shown  in  Fig.  lb.  From 
this  figure  we  can  obtain  the  following  useful  relations 


r  =  r2  sin  o 
dr 

r  =  r.  cos  p  . 
d « >  1 


(1) 

(2) 


The  following  notation  will  be  used  for  the  derivatives  with  respect  to  the 
independent  variables  0  and  <p 


&<>-<>’  <3> 

55 <  )  -  (  )'or(')  (4) 


2.  Stress  Resultants 

Integration  of  the  stress  components  through  the  thickness  h  yields  expres¬ 
sions  for  the  stress  resultants  acting  on  a  unit  length  of  the  section  0  =  const, 
or  9=  const,  (see  Fig.  2  for  notation  and  sign  conventions): 

h/2 

Np  =  /  +  z/r2)  dz  (5a) 

-h/2 

h/2 

Ne  =  /  ao<1  +  z/r^dz  (5b) 

-h/2 

h/2 

%  j  V1  +  z/r 2)dz  <5c> 

-h/2 

h/2 

NSp  =  /  +  z/rpdz  (5d) 

-h/2 
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(o)  STRESSES 


Fig.  2  Sign  Convention  for  Stresses  and  Stress  Resultants 
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h/2 

%  =  -  [  V  <1  +  z/r2)dz  (5c) 

-h/2 

h/2 

Q0=“/  T0z(1  +  z/rl)dz  (5f) 

-h/2 

h/2 

M0  =  -  j  ztTp  (1  +  z/r2)dz  (5g) 

-h/2 

h/2 

M0  =  -  J  zag(1  +  z/r^dz  (5h) 

-h/2 

h/2 

Mp0  =  "  /  ZT»0(1  "  z/r2,dz  (5,) 

-h/2 

h/2 

Mw  =  -  /  ZTot,  <*  +  z/ri» dz  (5J) 

-h/2 

The  present  analysis  is  for  a  thin  shell  (Assumption  n-c).  Therefore, 
h/r.  «  1  (i  =  1,2)  and,  in  Eqs.  (5),  we  neglect  the  terms  z/x^  in  comparison 
with  unity.  As  a  result,  =  N  )0  and  . 

3.  Constitutive  Equations 

The  Assumption  (Il-a)  that  the  shell  is  made  from  a  linearly  elastic  material 
enables  us  to  express  strains  in  terms  of  stresses  by  means  of  Hooke’s  law. 
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Since  normal  stresses  on  surfaces  parallel  to  the  middle  surface  are  neglected 
(Assumption  n-d),  Hooke's  law  takes  the  form 


%  "  E 

-  »a0) 

(6a) 

1  . 

€ti  E  '°0 

-  *"V 

(6b) 

■y  -  2(1  + 

r  pfl  E 

at 

<pe 

(6c) 

where  E  is  Young's  modulus  and  v  is  Poisson's  ratio. 

Solving  Eqs.  (6)  for  the  stresses  in  terms  of  the  strains,  we  find 


E 

c  = - ~(e  +  ve.) 

P  .  2 '  v>  O' 

1  -  V 

(?a) 

ae  =  ,  2  <fe  +  V 

1  -  V 

(7b) 

T  =  E  7 

>0  2(1  +  v)  \>0  * 

(7c) 

4.  Displacement  and  Rotation  Components 

The  circumferential,  meridional,  and  radial  displacement  components  of  a 
point  P  on  the  middle  surface  of  the  shell  are  denoted  by  u  ,  v  ,  and  w  , 
respectively  (see  Fig.  3a).  The  displacement  components  u  and  v  are  taken 
positive  in  the  direction  of  increasing  0  and  c>  ,  respectively,  w  is  positive 
when  it  points  away  from  the  center  of  curvature  of  the  meridian. 

In  the  present  analysis,  the  rotation  components  ,  and  will  be 

used  to  determine  the  contributions  of  prebuckling  quantities  to  equilibrium  of  a 
deformed  shell  element.  The  angle  of  rotation  of  the  normal  to  the  middle  sur¬ 
face  about  the  tangent  to  the  parallel  circle  at  P  is  denoted  by  (see  Fig.  3b). 
Now  because  of  the  Kirchhoff-Love  hypothesis  (Assumption  Il-d),  can  be 
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related  to  the  displacement  components  of  the  middle  surface.  To  establish  this 
relationship,  let  us  consider  the  displacement  of  two  neighboring  points  P  and 
Q  on  the  same  meridional  curve  of  length  r^h  (see  Fig.  4).  The  radial  dis¬ 
placement  of  Q  exceeds  that  of  P  by  an  amount  wd<  (see  Fig.  4a).  Conse- 
quently,  the  tangent  to  the  meridian  at  P  rotates  through  an  angle  wdo/r^dp 
and,  by  virtue  of  the  Kirchhoff-Love  assumption,  the  normal  at  P  rotates 
through  the  same  amount.  The  displacement  v  along  the  meridian  causes  the 
normal  at  P  to  rotate  through  an  angle  v/r^  (see  Fig.  4b).  The  circumferential 
displacement  u  does  not  contribute  to  the  rotation  of  the  normal  about  the  tangent 
to  the  parallel  circle  at  P  .  Hence,  by  adding  the  two  contributions  to  O  q  in 
accordance  with  the  sign  convention  shown  in  Fig.  3b,  we  obtain 


Proceeding  in  the  same  way,  we  can  show  that  the  angle  of  rotation  of  the 
normal  about  the  tangent  to  the  meridian  at  P  is  given  by 


*'  o 


^w*  -  u  sim 


(9) 


where  the  positive  sense  of  uu  is  shown  in  Fig.  3b. 

The  rotation  of  the  shell  element  around  the  normal,  denoted  by  u>  (see 
Fig.  3b),  is  not  clearly  defined  because  of  the  in-surface  shear  deformation. 

The  angle  of  rotation  ...  of  the  tangent  to  the  parallel  circle  around  the  normal 
differs  from  the  angle  of  rotation  uz2  °f  the  tangent  to  the  meridian  around  the 
normal.  Instead  of  using  an  average  rotation  u?  =  u;  .  +  w  ,  we  more 

Z  Z  i  Zm 
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properly  use  the  two  separate  rotations  and  jj ^  when  we  determine  the 
contributions  of  the  prebuckling  quantities  to  the  equilibrium  equations  of  a 
deformed  element.  We  let 


where  it  is  understood  that  o>  =  u>  .  if  the  prebuckling  quantity  under  consid- 
eration  acts  in  the  direction  of  the  tangent  to  the  parallel  circle,  and  that 
Cu'z  =  u  if  the  prebuckling  quantity  acts  in  the  direction  of  the  tangent  to  the 
meridian.  From  Ref.  13,  and  are  related  to  the  displacement 

components  of  the  middle  surface  by 


cos  $ 


5.  Equations  of  Kinematics 

Let  A  be  a  point  located  at  a  distance  z  from  the  middle  surface  and  let 
the  normal  through  A  intersect  the  middle  surface  at  the  point  P  .  From 
Eqs.  (VI-2)  Ref.  13,  the  strain-displacement  equations  for  point  A  are 


€4> 


va  :  WA 

’>(■  *  3 


(12a) 


U'a  +  VA  COS  0  +  WA  Sin  0 


'('  *  i) 


(12b) 


uA  cos  9  - 


’■(■  *  t)  '(‘  * « 


(12c) 
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f 


Once  again,  we  invoke  the  thin-shell  assumption  and  thus  neglect  z/r.  in  Eqs. 
(12).  Hence 


VA  +  WA 


uA  +  vA  cos  o  +  wA  sin  p 


UA  UA  C0S  9  "  VA 


pO  r. 


(13a) 


( 13b) 

(13c) 


The  right-hand  sides  of  Eqs.  (13)  contain  the  displacement  components  uA  , 
vA  ,  and  wA  of  point  A  .  Now,  the  Kirchhoff-Love  hypothesis  enables  us  to 
relate  uA  ,  vA  ,  and  wA  to  the  displacement  components  u  ,  v  ,  and  w  of 
point  P  on  the  middle  surface  and  to  the  rotation  components  and  . 
Since  and  can  be  expressed  in  terms  of  u  ,  v  ,  and  w  by  means  of 
Eqs.  (8)  and (9),  we  see  that  the  adoption  of  this  hypothesis  is  equivalent  to  re¬ 
ducing  the  problem  of  determining  the  displacements  of  the  shell  to  that  of 
determining  the  displacements  of  the  middle  surface.  The  Kirchhoff-Love 
hypothesis  implies  that  the  displacement  components  vary  linearly  through  the 
thickness  and  that  the  radial  displacement  component  wA  is  independent  of 
z  ,  i.  e. , 


UA  u  +  ~0Z  {14a) 

VA  v  v  "<)A  (14b) 


w.  -  w 
A 


(14c) 


( 


Hi 


By  inserting  into  Eqs.  (14)  the  expressions  for  and  given  in  Eqs.  (8)and  (9), 
we  obtain  the  relations  between  the  displacement  components  of  points  A  and  P  : 


=  u  +  (u  sinp  -  w')  ~  =  u(l  +  --  J-  w'~ 


(15a) 


vA  =  v  +  (v  -  w)  p-  =  v(l  + 


w.  =  w 


•  z 
~  WrT 


(15b) 

(15c) 


These  equations  are  the  same  as  those  derived  by  Flugge  (Ref.  13).  The  more 
general  elastic  law  derived  in  Ref.  13  was  simplified  for  a  thin  shell  by  re¬ 
placing  the  terms  (1  +  z/r j)  by  unity.  In  Ref.  13,  and  in  the  present  analysis, 
this  thin  shell  approximation  was  introduced  into  the  equations  defining  the  stress 
resultants  |Eqs.  (5)  J  and  into  the  strain-displacement  relations  |Eqs.  ( 12)  j . 
Reference  13  applied  the  thin  shell  approximation  to  Eqs.  (15)  and,  as  a  result, 
neglected  the  underlined  terms  in  Eqs.  (15).  For  example,  in  Eq.  (15b),  the 

term  v  —  was  neglected  in  comparison  with  v  .  However,  in  the  present 
rl 

analysis,  v/r ^  is  interpreted  as  being  a  part  v>f  the  rotation 


and,  as  such,  is  not  neglected.  Therefore,  the  underlined  terms  in  Eqs.  (15) 
are  not  neglected  in  the  present  analysis. 
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6.  Elastic  Law 


( 


From  Eqs.  (13)  and (14),  we  obtain  the  strains  ^  and  at  a  distance 

z  from  the  middle  surface  in  terms  of  the  displacements  u  ,  v  ,  and  w  at  the 
middle  surface  and  their  derivatives: 


*0  =  *0  -  2*0 

€e  =  €e  ‘  ZKe 

yq>e  =  y9v  =  ‘  z^2k00)  * 


where 

—  v  +  w 
■  rx 

-  u'  +  v  cos  0  +  w  sin  o 

(o  =  F 

—  U  U  COS  0  -  v' 

-  r1  r 

..  1  • 

<t>  =  "  r^ 

K0  =  “7^  -7^C0SP 
2k00  =  “  7^  ^0  "7  a'o  +  7  ^0  cos  0  » 


(16a) 

(16b) 

(16c) 


(17a) 

(17b) 

(17c) 

(17d) 

(17e) 

(17f) 


and,  from  Eqs.  (8) and (9), 


Ur 


w  -  v 


w*  -  u  sin© 


) 
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The  strains  ,  e0  ,  and  from  Eqs.  (16)  are  now  entered  on  the  right- 
hand  side  of  Eqs.  (7)  to  give  the  stresses  ,  a&  ,  and  in  terms  of  u  , 
v  ,  and  w  .  We  now  introduce  these  stresses  into  the  integrals  |Eqs.  (5)|, 
which  define  the  stress  resultants,  and  finally  obtain  the  elastic  law  for  a  shell 
of  revolution: 


N0  =  ***&  +  *««)  = 


No  - 


N 


00 

0 


0  ^  *Q‘ 


^0  +  *V 


pq..-..  v)y 

2  70O 


v  +  w  ,  /u'  +  v  cos*  +  w  sinp  \ 

=  dFT  ' - ' - '] 

"f 


u*  +  v  cos 0  +  w  sinp 


♦  *4^) 


DU  -  -) 

2 


U_  U  COS  0  -  V* 

r,  r 


=  K(  Kq  +  vkq)  =  K 


M0  =  K<V  "V 


Kl(-  Fw0  -  |  COS  0 


M00  =  K<x  -  l')* 


00 


(18a) 

(18b) 

(18c) 


(.±ue)+  -iu,ecoB» 

.  i  vj  t 

) +  "f  rfe  )] 

nj^Wl^UBiBej  +  l^sJjcosp 

■  -fcMl 

K(1  - 

# 

1  / w*  -  u  sinp  \  .  J_( W  -  v\ 

,2rA  r  /  2r\ri/ 


(18d) 


(18e) 


=  K(1  -  v) 


-  -l-(w’  -ru.8ln.l)  cos0 


2r 


(18f) 
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where  the  extensional  and  bending  stiffnesses  are  given  by 


D  = 


(19) 


and 


K 


Eh" 


12(1  -  v) 


(20) 


respectively. 

(Note  that  in  a  Donnell-type  analysis,  u  and  v  are  discarded  from  the 
formulas  for  k  ,  k  ,  and  k  \  . 

p  e  w 

The  elastic  law  given  by  Eqs.  (18)  is  the  same  as  that  derived  by 
Reissner  (Ref.  12)  and  Gravina  (Ref.  20)  who  used  the  methods  of  differential 
geometry. 

The  elastic  law  gives  six  equations  for  nine  unknowns  (N^  ,  N0  ,  , 

M  ,  ,  M  .  ,  u  ,  v  ,  w)  and  therefore  are  not  sufficient  to  determine  the 

p  V  0<7 

unknowns.  The  additional  equations  needed  are,  of  course,  found  from  the 
conditions  of  equilibrium  which  will  be  derived  in  the  next  chapter. 


i 


IV 


EQUATIONS  OF  EQUILIBRIUM 

In  deriving  the  equations  of  equilibrium  for  the  linear  bending  analysis  of 
shells  it  is  irrelevant  as  to  whether  the  forces  and  moments  are  assumed  to  act 
on  a  deformed  or  on  an  undeformed  element  of  the  shell.  In  shell  stability 
problems,  however,  the  equations  of  equilibrium  must  be  written  for  an  element 
of  the  deformed  (buckled)  shell.  These  equations  will  be  derived  for  a  general 
shell  of  revolution  in  this  chapter. 

1.  Equilibrium  Equations  for  the  Stability  of  a  Shell  of  Revolution 

In  this  section,  equations  of  equilibr  ium  are  derived  based  on  the  classical 
assumption  that  the  effects  of  prebuckling  rotations  may  be  neglected. 

Figure  1  shows  the  middle  surface  of  a  differential  element  of  the  shell  in 
the  prebuckled  state.  At  the  center  of  the  element,  point  P  ,  there  is  shown  an 
orthogonal  right-handed  system  of  axes  X  ,  Y  ,  and  Z  with  the  X  axis  in  the 
direction  of  the  tangent  to  the  meridian  at  P  ,  the  Y  axis  in  the  direction  of  the 
tangent  to  the  parallel  circle  at  P  ,  and  the  Z  axis  in  the  direction  of  the  outward 
normal  at  P  .  The  force  and  moment  equations  of  equilibrium  for  a  differential 
element  of  the  buckled  shell  will  be  written  with  respect  to  the  X  ,  Y  ,  and  Z 
directions.  The  terms  that  contribute  to  the  equilibrium  equations  may  be  divided 
into  four  groups  denoted  as  Group  0,  Group  1,  etc.  Group  0  contains  only  the 
prebuckling  quantities  (N0O  ,  ,  N0?)0  ,  p0  ,  p0  ,  pz)  .  Due  to 

prebuckling  equilibrium,  the  net  contribution  of  these  terms  to  the  equilibrium 
equations  for  the  buckled  shell  must  be  zero.  Group  1  contains  only  the  incre¬ 
mental  quantities  which  develop  as  the  shell  passes  from  the  prebuckled  state  to 
the  buckled  state,  h.  a  linear  stability  analysis,  only  terms  which  are  linear  in 
the  (infinitesimal)  incremental  quantities  are  retained.  These  terms  are  the 
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same  as  those  present  in  the  linear  bending  analysis  of  shells.  The  Group 
terms  arise  because,  in  the  buckled  state,  the  prebuckling  quantities 


^  act  on  a  deformed  element.  These  terms  con- 

P(p  *  p0  •  pz 

sist  of  products  of  a  (finite)  prebuckling  quantity  and  an  (infinitesimal)  incre¬ 
mental  displacement  or  its  derivative.  The  contributions  to  the  equilibrium 
equations  of  terms  which  are  linear  in  the  incremental  quantities,  i.e.,  the 
terms  in  Group  1,  2,  and  3,  are  given  in  the  next  three  subsections. 


In  the  determination  of  the  contributions  of  the  incremental  stress  resultants 


N0  ’  ‘  ’  %••••>  M04>  t0  the  equations  of  equilibrium  for  the  buckled  shell, 

it  is  irrelevant  as  to  whether  the  stress  resultants  are  assumed  to  act  on  a 


buckled  or  on  an  unbuckled  element.  This  is  because  the  equilibrium  equations 
for  the  buckled  element  differ  from  those  for  the  unbuckled  element  only  in  terms 
which  are  nonlinear  in  the  incremental  quantities,  and,  since  the  incremental 
quantities  are  infinitesimals,  these  nonlinear  terms  vanish.  Thus,  the  contri¬ 
butions  of  the  incremental  stress  resultants  may  be  obtained  from  a  consideration 
of  equilibrium  for  an  undeformed  element.  The  results  of  such  a  consideration 
follow  (Ref.  13): 


IFi  =  <rty  +  riN<te  -  rjN#  cos*.  -  rQp  (la) 

£F2  =  <rV'  +  rlNf>  +  rl%  C0S9  '  rl  Q0  sln  *  <lb> 

£F3  “  ‘  rlN0  sln  9  •  rNV  •  rl  %  -  <*%>’  <»<=> 

£Mi  =  <rMw>"  +  riMi  *  riM«o  cos  9  -  rri  %  <ld> 

sm2  =  -  <rM/  -  riM;0  +  riMe  cos  *  +  •  He) 


Note  that  a  common  factor  do  dp  has  been  omitted  in  Eqs.  (1).  This  factor 


will  also  be  omitted  in  the  sequel  when  other  contributions  to  the  equilibrium 
equations  are  obtained. 

1.2  Group  2:  Prebuckling  Stress  Resultants 

When  the  sh _1!  buckles,  the  prebuckling  stress  resultants  rotate  and  thereby 
develop  components  which  contribute  to  the  equilibrium  equations.  These  compo¬ 
nents,  which  consist  of  products  of  a  prebuckling  stress  resultant  and  an  incre¬ 
mental  rotation,  will  now  be  determined  for  each  of  the  prebuckling  stress 
resultants  N  ,  N0O  ,  ,  and  N/J0  . 


1.2.1  N  . 

pO 

When  the  shell  passes  from  the  prebuckled  state  to  the  buckled  state,  the 
meridional  force  N^rdfl  .  acting  at  point  of  the  section  p  =  pQ  (Fig.  1 
and  Fig.  2?),  participates  in  the  incremental  rotation  u,  and  therefore  develops 

U 

a  component 


K1  *  rN«,Ow©d0 


which  points  in  the  direction  of  the  outward  normal  (to  the  prebuckled  shell)  at 
.  Similarly,  the  meridional  force  N^rdO  +  (N^0rd9)*dp  ,  acting  at  point 
M2  of  the  section  p  =  pQ  +  dp  ,  develops  a  normal  (or  radial)  component 

Ki  ■  |Vd9  +  (fVord9>'H  K  *  “iH  =  rN0o“ed9  +  <rN0O<.0dO)'d0 

which  points  toward  the  center  of  curvature  of  the  meridian  at  point  Mg  .  The 
contributions  of  the  components  K  ^  and  K  *  to  the  equilibrium  equations  can 
easily  be  seen  from  Fig.  2b.  Thus,  due  to  the  difference  in  their  directions, 
Kj  and  K*  contribute  an  amount  to  the  equilibrium  of  forces  in  the 

{ 
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(c) 


Fig.  2  Contribution  of  N,- 


direction  of  the  tangent  to  the  meridian  at  P  ,  i.e. ,  in  the  X  direction;  and  due 
to  a  difference  in  the  magnitudes  of  and  K*  ,  there  is  a  contribution 
-<rN*o  l x)q )  to  the  equilibrium  of  forces  in  the  direction  of  the  normal  at  P  , 
i.e.,  in  the  Z  direction. 

Due  to  the  incremental  rotation  ^  0  ,the  meridional  force  N  ,.rd0  ,  at  the 

z2  pO 

section  p  =  pQ  ,  develops  a  component 


K2  rNp0^z2d0 


which  points  in  the  negative  Y  direction  (see  Fig.  2c),  and  the  meridional  force 

rN  Ado  +  (rN  do)  dp  ,  at  the  section  p  -  pA  +  dp  ,  develops  l  component 
01)  l) 


K2  ‘  |rV°  +  <rN,0d0)'dO|k2  +  -Z2d'"l  *  rN(,0-z2d°  +  <rNoO*Z2d°)'d'’ 


which  points  in  the  positive  Y  direction.  Since  both  K ^  and  K ^  act  in  the  Y 
direction,  there  is  a  contribution  +(rN  qu;^)*  to  equilibrium  of  forces  in  the 
Y  direction  as  a  result  in  the  difference  in  magnitudes  of  K ^  and  . 

The  stress  resultant  N  0  does  not  develop  any  components  due  to  the 

rotation  <*•  since  the  direction  of  N  .  is  parallel  to  that  of  u.  . 

p  pO  p 

Thus,  the  contributions  of  the  prebuckling  stress  resultant  f  to  the 
equations  of  equilibrium  of  forces  in  the  X  ,  Y  ,  and  Z  directions  are 

*rNpO^O’  +(rNP0  tcz2)  *  and  ~(rN00^)'  ’  respectively. 


1.2.2  N 


00 


The  hoop  force  N^r^dp  ,  acting  at  point  Cj  of  the  section  0  -  0 G  (Fig. 
3a),  participates  in  the  incremental  rotation  u.;i  and  develops  a  component 


S  “  rlN»oV' 


which  points  in  the  direction  of  the  outward  normal  at  .  Similarly,  the  hoop 
force  N^r^p  +  (N^r^pl’de  ,  acting  at  point  C2  of  the  section  0  =  oQ  +  dg  t 
develops  a  normal  component 

K3  '  hNeod*  *  <riNeod«'Hk+  Sd0| 

-  riN«o“e»d*>  +  (riN8oV,)d<’ 

which  points  toward  the  center  of  curvature  of  the  meridian  at  point  C2  .  From 
Fig.  3a  we  see  that  due  to  a  difference  in  their  directions,  the  forces  and 
K*  contribute  an  amount 


•<rlN90%8in,>) 

to  the  equilibrium  of  forces  in  the  Y  direction,  and  that  due  to  a  difference  in 
the  magnitudes  of  and  there  is  a  contribution 

-rl<N80“p>’ 

to  the  equilibrium,  of  forces  in  the  Z  direction. 

Due  to  the  incremental  rotation  u>zl  ,  the  hoop  force  N^r^do  ,  at  0  =  0Q  , 
develops  a  component 


K4  *  riNeo"zid* 

which  acts  along  the  tangent  to  the  meridian  at  0  =  0q  and  which  points  in  the 
direction  of  increasing  (Fig.  3b),  and  the  hoop  force  N^rjdp  +  (NgQrjdpl’d©  , 
at  the  section  0  =  0Q  +  d0  ,  develops  a  component 

K4  =  +  (rlNSOd^>'dHKl  +  "ild»l 

=  rlN80“zld8  +  <riN80“Zld5,)'d» 


which  acts  along  the  tangent  to  the  meridian  at  0  =  0Q  +  d0  and  which  points  in 

the  direction  of  decreasing  9  .  From  Fig.  3b  we  see  that  the  components  K4 

and  K*  contribute  an  amount 
4 

-ri(Neo“zi>' 


to  the  equilibrium  of  forces  in  the  X  direction  and  an  amount 


"riNeo“zi cos  ^ 


to  the  equilibrium  of  forces  in  the  Y  direction. 

The  stress  resultant  N0Q  does  not  develop  any  components  due  to  the  rota¬ 
tion  Ufy  . 

Thus,  the  contributions  of  the  prebuckling  stress  resultant  to  the 
equations  of  equilibrium  of  forces  with  respect  to  the  X  ,  Y  ,  and  Z  directions 

are  -riN„o  V"'*  -  r^^cos*  ,  awl  -r^N^)’  , 

respectively. 

1.2.3  Nrt  n  and  N 

0g>  0 _ flQO 

The  contributions  of  the  prebuckling  stress  resultants  N  ^  and  to  the 

equations  of  equilibrium  for  an  element  of  the  deformed  shell  can  be  obtained  in 
the  same  way  as  the  contributions  of  and  N0Q  .  Thus,  according  to  Figs. 

4  and  5,  we  obtain  for  equilibrium  of  forces  in  the  X  ,  Y  ,  and  Z  directions 
the  contributions 


-ritWz2C0Se  -  rNoe  o“o  -  <rNeoo“zi>  ’ 

"riNoaou'o  sin  9  *  ri(fWz2>'  •and 

-rl<N0(»0t“0 >'  ■  rlN0  00“z2sin0  •  <rtW'  +  rtWo"zl  ’  respectively. 


z 


Fig.  5  Contribution  of  ^ 
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1,2.4  Summary  of  Group  2  Contributions 

The  contributions  of  the  prebuckling  stress  resultants  Q  , 

and  to  the  equations  of  equilibrium  for  an  element  of  the  deformed  shell 
are 


SF?  -  -rNt,0U9  *  rl<NOo"zl>'  *  riNQo0OJz2  -  rN0flO“<, 
-(rN»eowzi>' 

ff2  -  *<rN0O“z2>’  *  rlNeOw0Bln«'  *  rlN90'Jzl  'OS0 
‘rlN90O“0  8ln*+  rl<N00O"z2)’ 


•(rN0Ou0)'  '  rl(N«OW0>'  ‘  ri(N00O‘V'  ‘  rlNe0O“z2  Sin  » 
■(rN00O“0)'  +  rN00o“zl  ’ 


1.3  Group  3:  Applied  Loads 

The  contributions  of  the  applied  loads  to  the  equilibrium  equations  will 
now  be  determined.  These  contributions  arise  because,  in  the  buckled  shell, 
the  applied  loads  act  on  an  element  which  has  been  deformed  by  the  incremental 
displacements. 

Let  the  components  of  the  applied  load  per  unit  area  of  the  middle  surface 
of  the  prebuckled  shell  be  denoted  by  p  ,  p.  ,  and  p  as  shown  in  Fig.  1. 

u  Z 

The  loads  p^  and  pfl  are  taken  positive  in  the  direction  of  increasing  q>  and 
$  ,  respectively;  p  is  positive  when  it  points  away  from  the  center  of  curvature 
of  the  meridian.  New  the  statical  approach  used  here  to  obtain  the  buckling 
loads  is  applicable  only  to  conservative  systems  (Ref.  21).  Hence,  the  load 
components  are  assumed  to  be  conservative;  for  example,  constant  directional 
loads  (dead  weight  loads)  or  hydrostatic  pressure  loads. 


Since  all  stability  problems  involve  considerations  of  the  deformed  struc¬ 
ture,  it  is  necessary  to  specify  precisely  the  character  of  the  applied  loads;  i.e. , 
the  way  the  applied  loads  behave  as  the  shell  deforms.  Thus,  for  example,  the 
stability  equations  for  a  constant  directional  pressure  loading  will  differ  from 
those  for  a  hydrostatic  pressure  loading.  It  is  conceivable  that  small  changes 
in  the  nature  of  the  applied  loads  might  have  an  appreciable  effect  on  the  magni¬ 
tude  of  the  buckling  load;  indeed,  for  buckling  of  a  ring  the  difference  amounts 
to  33%  (Ref.  22).  Two  types  of  pressure  loadings  will  be  considered  here:  con¬ 
stant  directional  pressure  loading  (dead  weight  loading)  and  hydrostatic  pressure 
loading.  A  constant  directional  pressure  loading  is  such  that  the  total  force 
acting  on  a  shell  element  does  not  vary  in  magnitude  or  in  direction  as  the  ele¬ 
ment  deforms.  Thus,  a  constant  directional  force  does  not  contribute  (explicitly) 
to  the  equilibrium  equations.  A  hydrostatic  pressure  loading  is  such  that  the 
total  force  acting  on  an  element  of  the  shell  is  always  proportional  to  the  actual 
size  of  the  element  and  is  always  directed  normal  to  the  element.  Then,  due  to 
the  incremental  rotations  and  ,  the  hydrostatic  pressure  force  rr ^pzd0d0 
acting  on  an  element  of  the  shell  develops  the  components  rr^u^dfldp  and 
rr1Pzaipd0do  which  point  in  the  positive  X  and  Y  directions,  respectively.  In 
addition,  due  to  the  stretching  of  the  middle  surface  during  buckling,  there  is  a 
component,  +rr.p  (7  +  7  )d0d0  which  points  in  the  positive  Z  direction. 

Thus,  the  explicit  contributions  of  the  applied  loads  to  the  equations  of 
equilibrium  for  an  element  of  the  deformed  shell  may  be  written  as 

SFi  =  +Vrri‘V'9) 

IF2  =  *5ph(rrlPz-0> 

IF3  =  +6ph  rr!Pz<70  +  V 


(3a) 

(3b) 


(3c) 


where 


1 ,  for  hydrostatic  pressure  loading 
^ph  0,  for  constant  directional  pressure 


(4) 


1 . 4  Summary  of  Results 

By  adding  the  contributions  given  by  Eqs.  (1)  through  (3),  we  arrive  at  the 
following  equations  of  equlibrium  for  an  element  of  the  buckled  shell: 


2F,  = 


<rV  +  riN<i„  -  riNe  008  p  +  |-rS,  -  rN»o“e  '  rN„„o%l 

+  |"ri(Neouzi>'  -  <rIU“zi>'  -  riNS(, 0^7.2  008  pl 
+  sphrrlpz^  =  0  <5a> 

<rV  +  riN«  +  ri%  cosp  +  l‘riQo  3ino  -  riN«0Oa'csin'( 

'riWH  +  l"riNeo“zi  co8p  +  (rNeo“z2>' 

+  rl<Ne(,0“z2>']  +5phrripza«  =  0  (5b) 

2F3  =  -r,H  sinp  -  rN0  -  -  <*}/  -  (rN^)'  - 


TF„  = 


-rl(NOOu,P)'  ‘  (rNPCoV’  +  Koo^’zl  *  rlNe W01z2  8in  8 


+  6ph  rrlpz‘v  +  V  =  0 


EM,  = 


ZMn  = 


(rM  J  +  r,M'  +  r,M  cos©  -  rr,Q  =0 
P0  1  9  1  OP  If) 

-(rM0)  -  riMoo+  T^ecos(t>  +  rrxQ  =  0 


(5c) 

(5d) 

(5e) 


If  effects  of  rotations  around  the  normal  are  negligible,  then  the  terms  in  the 
straight  brackets  in  Eqs.  (5)  are  omitted.  In  a  Donnell  type  analysis,  the  terms 
in  the  braces  and  brackets  in  Eqs.  (5)  are  neglected. 


2.  Nonlinear  Equations  of  Equilibrium  for  a  Shell  of  Revolution 


The  stability  equations  presented  in  the  preceding  section  were  based  on 
the  classical  assumption  that  the  prebuckling  rotations  could  be  neglected.  In 
this  section,  we  present  stability  equations  which  include  the  effects  of  prebuckling 
rotations.  These  equations  are  obtained  through  specialization  of  the  nonlinear 
equations  of  equilibrium  which  will  be  derived  here.  To  obtain  the  stability 
equations ,  we  replace  each  unknown  quantity  O  In  the  nonlinear  equations  of 
equilibrium  by  (  )Q  +  (  )  ,  where  (  )Q  represents  the  prebuckling  value  of  the 
quantity  and  (  )  represents  the  incremental  value  of  the  quantity  which  develops 
during  buckling.  Then  in  the  resulting  equations,  the  terms  containing  only 
prebuckling  quantities  may  be  subtracted  out  by  virtue  of  prebuckling  equilibrium, 
and  nonlinear  terms  in  the  (infinitesimal)  incremental  quantities  may  be  dis¬ 
regarded.  The  result  is,  of  course,  the  linear  stability  equations.  Actually, 
only  the  nonlinear  force  equilibrium  equations  will  be  presented  here.  The 
nonlinear  moment  equilibrium  equations  can  be  obtained  from  considerations 
identical  to  those  used  in  arriving  at  the  nonlinear  force  equilibrium  equations . 
However,  for  the  case  of  a  thin  shell,  the  nonlinear  terms  in  the  moment  equi¬ 
librium  equations  are  negligible. 

The  stress  resultants,  rotations,  and  strains  which  develop  as  the  shell 
deforms  due  to  the  applied  loads  p .  ,  pu  ,  and  p  are  denoted  by  N  , 

<p  v  Z  (p 

. % . S60  "zl  ’  "z2  •  ’  “d  h  •  These 

quantities  are  zero  when  the  shell  is  unloaded.  Linearized  relations  between 

these  quantities  and  the  displacement  components  of  the  middle  surface  can 
be  obtained  from  Eqs.  (III-8),  (III -9),  (HI- 11),  and  (IH-18)  providing  the  dis¬ 
placement  components  appearing  in  these  equations  are  measured  from  the 
unloaded  shell. 

Figure  6  shows  the  middle  surface  of  a  differential  element  of  the  unloaded 
sheli.  At  the  center  of  the  element,  point  P*  ,  there  is  shown  an  orthogonal 


right-handed  system  of  axes  X  ,  Y  ,  and  Z  with  the  X  axis  in  the  direction  of 
the  tangent  to  the  meridian  at  P*  ,  the  Y  axis  in  the  direction  of  the  tangent  to 
the  parallel  circle  at  P*  ,  and  the  Z  axis  in  the  direction  of  the  outward  normal 
at  P*  .  The  equations  of  equilibrium  for  a  differential  element  of  the  deformed 
shell  will  be  written  with  respect  to  the  X  ,  Y  ,  and  Z  axes.  Let  M*  denote 
the  point  of  intersection  of  the  meridian  through  P*  and  the  upper  parallel 
circle  (see  Fig.  6).  At  the  point  M*  let  the  directions  of  the  tangent  to  the 
meridian,  the  tangent  to  the  parallel  circle,  and  the  normal  be  given  by  the  X*  , 
Y*  ,  and  Z*  axes,  respectively.  Now  due  to  the  applied  loads,  the  line  M*A  , 
which  is  tangent  to  the  meridian  of  the  unloaded  shell  at  M*  ,  acquires  the  new 
direction  X**  as  shown  in  Fig.  7a.  The  meridional  force  rN^d0  ,  at  point 
M*  of  the  section  c  -  C.  ,  is  defined  to  act  in  the  X**  direction.  Then  ac- 

i  V 

cording  to  Fig.  7a,  the  components  of  N^  along  the  X*  ,  Y*  ,  and  Z*  direc¬ 
tions  are 


rN.  cos  ul  cos  u,'  d0  , 

V  u 

rNp  cos  -q  sin  uJz2  d0  ,  and 

rN0  sin  cl’gd0  ,  respectively.  These  components 
are  shown  in  Fig.  7b.  The  components  of  the  meridional  force  rN^d0 
+  (rN^d0)  d$  ,  which  acts  at  point  M*  of  the  section  *  =  +  dp  ,  are  also 

shown  in  Fig.  7b.  Next,  with  sin  Z.  — •  Z.  and  cos  Z.  -*  1  ,  and  from 
Fig.  7b,  the  contributions  of  the  stress  resultant  $0  to  the  equilibrium  equations 
written  with  respect  to  the  X  ,  Y  ,  and  Z  directions  are  (rN^)*  -  rN^  u30  , 

—  •  — r-  « 

(rN0£z2)  ,  and  -rN^  -  (rN0Gfl)  ,  respectively. 

The  contributions  of  the  other  stress  resultants  may  be  obtained  in  the 
same  way  as  the  contributions  of  were  obtained.  Thus,  according  to 
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Figs.  8  to  10,  the  nonlinear  equations  of  equilibrium  for  a  differential  element 
of  the  deformed  shell  are 


ZF 


1 


ZF 


3 


sm2 


<rV  +  (rlV'  ‘  rlV08  0  -  rQ?  -tvfue- 

'  (rR09C2i>'  -  <riITeszi)'  •  rlfr90Cz2  008  * 

+  rrlp®  *  6ph(rrlpz“e)  =  0 

<rlV  +  <rR0O»’  +  rl57W  cos  *-  r1^9  sin  *  *  rlT,W°98to  * 


(6a) 


-  rlF0  °0Sln  *  +  <rlS90°z2)'  +  <rVz2»' 

‘  rlK0Czl  008  *  *  rrlp0  +  VVzV  =  0  <6b> 

-rS»  -  ri% 8inp  -  (rV  -  <riV  -  <rfyv  '  iTlW 

-  <rlJ,#V'  '  <riV*>’  +  r%e°zi  -  rifVz28inp 

+  rrlpz  +  6phrrlpz(l«  +  V  =  0  (6c) 

(rM00)  +  (r^)’  +  cos  p  -  rr^  =  0  (6d) 

*(rB0)‘  -  (r^)'  +  rjM0  cos  <p  +  rr^Q^  =  0  .  (6e) 


The  equations  governing  the  stability  of  a  shell  of  revolution  may  now  be 
derived  from  the  nonlinear  equations  of  equilibrium.  In  Eqs.  (6),  we  let 


N0  N<*>0  +  N0  ’  *  *  ‘  ’  M00  M00O  +  M00  * 


=  w  +  , 

H>  VO 


(7a) 


(7b) 


where,  for  example,  is  a  prebuckling  quantity  and  is  the  increment 

in  this  quantity  which  develops  during  buckling.  After  insertion  of  Eqs.  (7)  into 
Eqs.  (6),  the  terms  containing  only  prebuckling  quantities  may  be  subtracted 
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Fig.  8  Components  of  N 


(o) 
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out  by  virtue  of  prebuckling  equilibrium,  and  nonlinear  terms  in  the  incremental 
quantities  may  be  neglected  since  the  incremental  quantities  are  considered  to 
be  infinitesimal.  This  procedure  results  in  the  replacement  of  the  nonlinear 
equations  of  equilibrium  by  linear  equations  of  stability.  These  equations  are 


<rV'  +  <riV'  -  riNo  008  0  -  rQs  -  r<N„0‘V  vw 

r*N0co“e  *  N00ti»o)  '  *rN®0o“zi  *  rN0ou'zio’ 

'  *riN0o“zi  +  riNe“zio* 

-  rl<N00OUz2  *  Nm“z20>  C0S  0  +  rrl5phPz‘1)  =  0  <8?> 

<riV  +  <rV’  +  riNe0cos«>  -  ri V,n» 

■  ri(,W'0  *  N»»<W  sln0  -  ri<Noo“*  +  Np<W  sin0 
+  <rlN«0(>“z2  +  riN0o“z2O)  +  (rNpO<1,z2  +  rN0U’z2O) 

'  rl(N«0^zl  +  VzlO>  C0S  0  +  rrl  VzS  r  0  <8b> 

*rV  •  rlN0  sin0  -  (rV  '  <rlQ0>'  '  <rN»0“0  +  rV'eO>' 

-  <rlN00Oc-0  +  rlN00IV'  -(rN00O%  *  rN00u'«o' 

-  <rlN0O%  +  rlVpO>'  +  r(No0Oazl  *  N00“zlO» 

■  rl(N00O“z2  +  N0 0Wz2O>  Sin0  *  rrl4phpz  <f0  +  r0>  "  0  <8c> 

(rM  )  +  r  M ’  +  r.M  cos  p  -  rr.Q  =  0  (8d) 

1  0  1  Op  it) 

-(rM^)  -  rMl^  +  r^M^  cos  <t*  +  rr^Q^  ~  0  •  (8c) 


When  prebuckling  rotations  are  neglected,  these  equations  reduce  to  the 
previously  derived  equations  of  stability  lEqs.  (5)] . 

Nonlinear  equations  of  equilibrium  for  a  general  shell  have  been  derived  by 
Sanders  (Ref.  23)  and  Kempner  (Ref.  24).  In  Ref.  23,  the  derivation  was  carried 
out  in  tensor  form  whereas  Ref.  24  used  a  variational  approach.  The  results 
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of  the  present  analysis,  Eqs.  (6),  agree  with  those  of  Refs.  23  and  24  except  for 
terms  which  involve  rotations  around  the  normal.  Kempner  did  not  include  such 
terms  whereas  Sanders  used  an  average  rotation  around  the  normal  given  by 
u.1  =  u:  +  u;  The  present  analysis  uses  u>.  and  w  9  for  t1^  rotations 
about  the  normal  of  the  tangent  to  a  parallel  circle  and  the  tangent  to  the  meridian. 
Thus,  the  discrepancy  between  the  present  results  and  those  of  Sanders  is  due, 
in  part,  to  the  use  of  different  expressions  for  the  rotation  around  the  normal. 

3.  Equations  of  Equilibrium  for  a  Cylinder 

The  equilibrium  equations  (5)  will  now  be  specialized  for  a  circular 
cylindrical  shell  subjected  to  a  uniform  external  pressure  p  ,  and  an  axial 
compression  at  the  edges.  The  axial  force  per  unit  length  of  circumference 
is  denoted  by  P  .  The  sign  conventions  for  the  coordinates  (x,  d)  ,  dis¬ 
placements  (u,  v,  w)  ,  and  stress  resultants  are  shown  in  Fig.  11  (compare 

with  Figs.  Ill— 2  and  III-3).  Thus,  in  the  equilibrium  equations,  we  let  <p  —  ^  , 

“3  w 

r  d<P  dx  ,  r  «  ,  r  -♦  a  ,  and  r  *  a  .  Also,  we  let  u>  =  - — 

1  1  6  a  0X  | 

1  V  m  _  Z 1  <L*L  m  =  —  C  -  —  +  1  9v  C  - 

a  'as  "  V)  ’  zl  a  30’  z2  9x  »  0  a  a  90  ’  dx  ’ 

N^0  -  -P  ,  NQ0  =  -pa  ,  N0qo  =  0  *  and  5ph  =  1  (a88Ume  P  is  a  hydrostatic 
pressure).  This  yields  a  set  of  equilibrium  equations  which  are  identical  to 
those  derived  by  Fliigge  in  1932  (Ref.  2): 


JL  n  *  J.  JLn 

8x  *  x  a  dd  dx 


aw  i  a  u 
8X  ’  3  302 


(19a) 


—  -2-M.  _  JL  m 
a  dd  b  ax  x0 


-  Pa 


(19b) 


-  ~r  M  +  M  n 

a  „„2  d  9x9 0  xO 

dd 


-5-  M  +  N 
£  x  0 


_u_  +  w  +  _i  a2w) 

)x  +  a  a  ’  as2/ 


+  Pa 


(19c) 
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(a)  COORDINATES  8  DISPLACEMENTS 


(b)  STRESS  RESULTANTS 


Fig.  11  Conventions  for  Coordinates,  Displacements,  and  Stress 
Resultants  of  a  Cylinder 
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4.  Equations  of  Equilibrium  for  a  Sphere 


The  equations  of  equilibrium  for  axially  symmetric  buckling  oi  a  sphere 
loaded  by  a  uniform  hydrostatic  pressure  p  are  obtained  from  Eqs.  (5)  by 
lecting  =  a  (radius  of  sphere),  r  =  a  sin  p  ,  =  ~  -  v)  , 

*o  =  STm?*''008  *  +  wsln*>  •  '0  =  7<w  +  v‘>  •  N„o  =  -  of  -  Neo  = 

and  )  =  u  -  N =  wzi  =  w Z2  ~  ~  0  ‘  yields  a  set  °f  equilib¬ 

rium  equations  which  are  identical  to  those  derived  by  Fliigge  (Ref.  13): 

(N^  sin  p)  -  N0  cos  p  -  sin  p  +  -^(w-  "  v)  sinp  =  5  (10a) 

(N^  +  Ng)  sin  q  +  (Q^  sin  p)  +  -|-(v  sin  p  +  v  cos  p 

•  . . 

+  4  w  sin  p  +  w  cos  p  +  w  sin  p)  =0  (10b) 

(Mp  sin  p)  -  M0cos  p  -  a  Q^sinp  =  0  .  (10c) 
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SHELLS  OF  REVOLUTION  UNDER  AXIALLY  SYMMETRIC  LOADS 


In  this  chapter,  and  in  the  sequel,  only  axisymmetrically  loaded  shells  of 
revolution  are  considered.  Also,  the  eifects  of  prebuckling  rotations  are  neglected. 
Since  the  loading  is  axially  symmetric,  the  space  variables  in  the  partial  dif¬ 
ferential  equations  governing  the  stability  of  a  shell  of  revolution  can  be  sepa¬ 
rated.  After  the  separation  of  variables,  three  equations  in  terms  of  the  three 
displacement  components  are  obtained  through  combination  of  the  equilibrium 
equations  and  the  elastic  law. 

1.  Separation  of  Variables 

The  stability  equations  for  a  shell  of  revolution  are  given  by  a  system  of 
partial  differential  equations  with  variable  coefficients.  For  the  case  of  axially 
symmetric  loading,  the  coefficients  in  these  equations  are  independent  of  the 
circumferential  coordinate  6  (see  Fig.  HI-1  for  notation).  Consequently,  it  is 
possible  to  separate  variables  and  thus  replace  the  system  of  partial  differential 
equations  by  a  system  of  ordinary  differential  equations .  Such  a  separation  of 
variables  is  effected  by  means  v,f  the  following  Fourier  series  representation 
for  the  incremental  quantities: 


.  £ 


V  =  2  vn<0)cos  ne 


00 

u  =  ^  un  (0)  sin 


w  -  J  wn  (0)  cos 


%  ~  '  2  *W0)8in 


N0  =  2_N0n<0)cos  nd 


=  %  =  2  M00n  8in  n0 


w 

Ne  “  2  N0n  <♦>  cos  ne 


Qe  -  2  Q0n  (0)  stn  ne 


w 

=  ]>  M0n<0)cos  nfl 


oo 

=  2  °*n (0)  8in  n0 


!0  =  J  M0n^cos  n0 


“zl  =  2  "zln  <*)  sl” 


S»  =  2  Q*n<0,cos 


“z2  =  Z  “z2n  <0*  sin  "9 


"e  '  2  “«n  <*> 008  ■*» 


°  2  %n  <(,)  cos 
n=0 


£s  ‘  2  €<m  <0> 008 
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These  expressions  represent  the  general  solution  of  the  stability  equations  for 
the  cases  in  which 

(i)  the  shell  is  complete  in  the  circumferential  direction,  i.e. ,  0<  0  <  2 x 

(ii)  the  boundary  conditions  on  the  edges  0  =  ,  —  for  a  partial 

toroidal  shell  are 

w  =  =  v  =  N0  =  0  .  (2) 


For  the  first  case,  the  requirement  that  the  incremental  quantities  should  be 
periodic  functions  of  0  is  satisfied  by  each  term  of  the  Fourier  series  in  Eqs.  (1). 
For  the  second  case,  the  boundary  conditions  ^Eq.  (2)j  are  satisfied  by  each 
term  of  the  Fourier  series  obtained  by  assigning  values  to  n  in  Eqs.  (1)  as  follows: 


n  =  n* -g-  ,  n*  =  1,  3,  5,  . . . 


(3) 


2.  Equilibrium  Equations 

Insertion  of  the  Fourier  series  representations  of  the  incremental  quantities 
| Eqs.  (l)j  into  the  equations  of  equilibrium  |Eqs.  (IV-5)j  results  in  the  following 
equations  of  equilibrium  for  each  value  of  n  : 


ZF 


2F 


in  =  <rV  *  nrlNWn  -  rlNon  cos*  -  T%n  '  rN„Owen 
-  nrlN0Owzln  +  rrlS  ph  pzw0n  ’  0 
2n  '  <rtW'  •  nrlN»n  +  sin  O 

-  rlN»%  6ln  ?  -  rlNw“zlncosp 

4  <riN0O  008  *  4  rV  "z2n  4  <rV  “z2n 

4  rrl*  phPz^n  =  0 


(4a) 


(4b) 
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SF3n  =  rV  *  rlNen  8in  *  +  “rlQ<m  +  (rV 

+  <riN0O  ccs0  +  rV"«n  +  rN»owm 

+  nrlNfl0Vrrl*phPz(%n+  V 

=  0 

(4c) 

JMln  =  <r*W‘  *  "An  +  rJ*W!0a  *  *  "lSjn 

=  0 

<4d) 

IM2n  =  <r*V‘  +  "lM*»n  ‘  r  An  008  ♦  *  rrl%n 

=  0 

(4e) 

Elastic  Law 

The  following  relations,  which  are  obtained  from  the  elastic  law 

[Eqs.  (in-18)], 

can  be  used  to  express  the  equilibrium  equations  in  terms  of  the  displacement 

components  u  ,  v  ,  and  w 
n  ’  n  ’  n 


5<rV 

£><rlN0n> 
D  *rlN0$n* 
K  ^rlM0n^ 

K  *rlM0$nJ 


+ 


v  + 
n 


(1  -  v)  n 


w  + 


n 


~U  -  v)  n  - 


(5c) 
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where,  for  brevity,  we  have  introduced  the  notation: 


Xj  =  sin  P  ;  yj  =  cos0  ;  and  y2  =  cos  2b  (6) 


4.  Rotation  and  Strain  Components 

By  substitution  of  Eqs.  (1)  into  Eqs.  (in-8),  (III— 9) ,  (1U-11),  and  (III-17),  we 
obtain  the  following  expressions  for  the  rotation  and  strain  components: 


u> 


0n 


c o 

On 


+ 


u 

n 


CO 


zln 


ii 

r 


v 

n 


co 


z2n 


n 


(7a) 


(7b) 

(7c) 

<7d) 


(7e) 


(7f) 


5.  Stability  Equations  for  the  Axially  Symmetric  Loaded  Shell  of  Revolution 

For  a  thin  shell,  the  effect  of  the  transverse  shear  forces  Q.  and  QA  in 

<pn  0n 

the  first  two  equations  of  equilibrium  (Eqs.  (4)]  may  be  omitted.  (Also  we  note 
that  Steele  (Ref.  25)  has  pointed  out  that  the  retention  of  these  terms  in  the  first 
two  equations  of  equilibrium  is  inconsistent  with  the  use  of  an  uncoupled  elastic 
law.  ]  The  moment  equations  of  equilibrium  (Eqs.  (4d-e)l  may  be  used  to  elimi¬ 
nate  Q^n  and  from  Eq.  (4c),  and  the  stress  resultants,  rotations,  and 


51 


strains  may  be  eliminated  from  the  equilibrium  equations  through  insertion  of 
Eqs.  (5  and  7)  into  Eqs.  (4).  In  this  way,  we  finally  obtain  a  set  of  three  ordinary 
differential  equations  for  the  three  displacement  components  u^p)  ,  vn($)  ,  and 

wnM  : 

<hl  +  °1>  un  +  <h2  +  c2>  “n  +  <h3  +  °3>  Un  +  (h4  +  c4>  vn  +  <hS  +  c5>  % 

+  <b6  *  c6>  Vn  +  (h7  +  c7>  Vn  +  (h8  +  c8>  wn  +  (h9  +  V  wn 

(8a) 

(fl  +  al>  un  +  <f2  +  a2>  un  +  <*3  +  a3>  “n  +  (f4  +  a4>  vn  +  (f5  +  a5>  vn 
+  <f6  +  V  \  *  <f7  +  V  \  +  (f8  *  a8>  Wn  +  <f9  +  V  "n 

+  <f10  +  a10>  %  +  <fll  +  all>  "n  +  <f12  +  a12>  "n  *  0  <8b> 

<gl  *  bj)  un  ♦  (g2  ♦  b2)  un  +  <g3  +  b3)  Un  ♦  <g4  +  b4)  vn  +  <g5  -  b5)  vn 

••  •••  • 

+  <g6  +  b6>  Vn  '  <g7  +  b7,vn  *  <g8  +  b8>  wn  +  (g9  +  V  Wn 

+  <g10  *  b10)wn  +  <gll  +  VVn  +  (g12  +  =°  <8c> 


+  <h10  +  c10>  Wn  +  (hU  +  Cll>  w.,  +  <h12  +  C12>  wn  =  0 


where  the  coefficients  a,b,c,ftg,h,  which  depend  on  ,  ,  are  given 
by  the  following  formulas: 
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For  convenience,  a  subscript  n  was  omitted  in  the  expressions  for  the  coeffi¬ 
cients  a  ,  b  ,  c  ,  f  ,  g  ,  and  h  in  Eqs.  (8)  and  (9).  We  note  that  the  functions 
a (<p)  ,  b(<t>)  ,  and  c(0)  in  Eqs.  (8)  depenu  on  the  prebuckling  quantities  , 
NT0O  ,  and  py  whereas  f(6)  ,  g(<fc)  ,  h(0)  do  not. 
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VI 


TOROIDAL  SHELL  UNDER  EXTERNAL 
PRESSURE -THEORETICAL  ANALYSIS 


1 .  Stability  Equations  for  a  Toroidal  Shell  Under  External  Pressure 

As  an  application  of  the  theory  developed  in  the  preced  \g  chapters,  the 
equations  governing  the  stability  of  a  general  shell  of  revolution  will  now  be 
specialized  for  a  toroidal  shell  subject  to  a  uniform  hydrostatic  pressure  p  . 

The  notation  for  a  toroidal  shell  with  a  circular  meridian  is  shown  in  Fig.  1. 
The  radius  of  curvature  of  the  meridian  r^  is  denoted  by  a  ,  and  the  distance 
between  the  center  of  the  circular  cross  section  and  the  axis  of  revolution  is 
denoted  by  b  .  Note  that  the  shell  geometry  and  the  applied  loading  are  sym¬ 
metric  about  the  plane  A-A  (see  Fig.  1).  Hence,  it  is  expected  that  the  buckling 
pattern  will  be  either  syirmetric  or  antimetric  about  this  plane.  Accordingly, 
it  is  convenient  to  use  instead  of  <p  the  coordinate  ip  measured  from  the 
plane  A-A  as  shown  in  Fig.  1.  The  independent  variable  ip  is  related  to  the 
colatitude  (p  used  in  Chapter  V  by 

*  =  *  +  y  .  (1) 

From  Eqs.  (1)  and  (III-4)  we  find 


-7  ( 

)  = 

7-  (  ) 

=  (  )'  • 

(2) 

d<p 

For  brevity,  we  introduce 

the  notations 

S  = 
m 

;  sin 

ml  ,  m 

=  1.2,.. 

(3a) 

C  = 

m 

:  cos 

mi/  ,  m 

=  0,1,2,... 

(3b) 

m 


Fig.  l  Notation  for  a  Toroidal  Shell 


59 


Therefore,  the  notations  sin  (b  ,  cos  <b  ,  and  y ^  -  cos  2<£ 

used  in  the  preceding  chapter  are  replaced  by 


X1  "  C1  =  cos  ^ 

(4a) 

yl  =  ”S1  =  ~sin  ^ 

(4b) 

y2  =  "C2  K  "COS 

(4c) 

From  Fig.  1,  we  find  that 

r  -  b  +  a  cos  4  , 

which  can  be  rewritten  as 

X  „  a  *  C 
a 

(5) 

where  0  is  a  nondemens ional  geometric  parameter 

a  -± 

a 

(6) 

The  prebuckling  stress  resultants  and  N  may  be  obtained  from 


Ref.  (13): 


J. A  -  -  v — — : — -  (b  +  -|-sin  0) 
bO  b  +  a  sin  0  2 


N 


00 


JE2. 

2 


We  introduce  a  nondimensional  load  parameter 


(7) 


GO 


where  E  is  Young's  modulus  and  h  is  the  shell  thickness.  Then  the  prebuckling 
stress  resultants  may  be  written  as 


DN«0  =  -i'1  -  >'2)f(2a+  C1,X  (8a) 

TjN80  =  '  (8b) 


where  D  is  the  extensional  stiffness  of  the  shell  jsee  Eq.  (m-19)j. 

We  may  now  obtain  the  stability  equations  for  a  toroidal  shell  through 
specialization  of  the  equations  for  a  general  shell  of  revolution  |  Eqs.  (V-8)J. 


(iJ  •  E«- 


(V-8b)  by  -7  ,  and  Eq.  (V-8c)  by 

a 


We  multiply  Eq.  (V-8a)  by 
/r\2 

(— j  in  order  to  remove  any  dependence  on  t  in  the  denominators  of  the  coef¬ 
ficients  in  the  stability  equations.  Next,  we  divide  the  resulting  equations  by 
D  and  introduce  the  nondimens ional  geometric  parameter 


.  K  _  1  /h\ 

k  ^2  12\a/  *  <9> 

Finally,  we  substitute  los.  (i)  through  (9)  into  Eqs.  (V-8)  and  (V-9)  and  obtain 
the  following  stability  equations  for  a  toroidal  shell  subject  to  a  uniform 
external  pressure: 

(h1  +  c,  X )  uQ  ( ip )  +  <h2  +  c2  X)  u^  (ip)  +  (h3  +  c3  X  )  un  (ip) 

+  (h4  +  C4A)  vn  (ip)  +  (hg  +  cg  *)  vn  (<P) 

+  +  c6*^n^  4  (1*7  +  C7X> 

+  (h8  +  V)wn(W  +  (h9  +  °9X)  VW 

+  (h10  +  C10X)  VW  +  (hll  *  C11X>  V*> 

+  (^12  +  ci2X)w*n(^)  =  0  (10a) 
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(fx  +  axx)  un  (ip)  +  <f2  +  a2X)  un  (ip)  +  (f3  +  a3X)  iin  (if>)  +  (f4  +  a4X)  vn(if) 

•  t •  • 

+  (f5  +  agx)  vn  (if>)  +  (f6  +  agx)  vn  (if,)  +  (f?  +  a?X)  vn  fy) 

+  <f8  *■  a8X)  wn  m  +  <f9  +  V>  "n  <*>  +  (f10  +  a10X)  "n  <*> 

•••  •••• 

+  <fH  +  anx>  wn  W  +  <fi2  +  a12X)  wn  W  =  0  <10b) 

(gx  +  bjX)  un  (ifi)  +  (g2  +  b2X)  un  (ip)  +  (g3  +  bgX)  u'n  (tf) 

+  (g4  +  b4X)  vn  (ip)  +  (g5  +  bgX)  vn  ty)  +  (gg  +  bgX)  vn  (if>) 

•••  . 

+  (g?  +  b?X)  vn  (ip)  +  (gg  +  bgX)  wn  (ip)  +  (g9  +  bgX)  wn  (ip) 

•  •  •  •  •  ••  •• 

+  <K10  +  b10X)  wn  <*>  +  <«11  +  bnX>  wn  m  *  (g12  +  b12X)  wn<*>  =  0 

(10c) 

In  Eqs.  (10)  the  coefficients  a  ,  b  ,  c  ,  f  ,  g  ,  and  h  ,  which  are  now  different 
from  those  defined  in  Eqs.  (V-9),are  given  by 

hx  =  nk|  n2Cx  -  4C^Sy  -  (a+  CJ  C2  -  (2  -v)(a+  Cp  C2 
+  (1  +v)(a+  CpSj2  +  (a+  Cj)2  cj 

+  n|  (a  +  C])2C1  +  v<a  +  C^3]  (11) 

h2  =  nk[-3(a+  CpCjSj  +  2(c+  C^2  Sj  ] 
h3  =  -nk  («  «■  Cj)2  Cx 

h4  =  k[  +(a  +  C.JSj3  -  n2(a+  CjJSj  +  2(a+  Cj)2 

-  „(a  +  Cj)3  sj  -  (a  +  C^2  CjSj  -  v(o  +  C{)3  Sj 
h5  -  k  Jn2  (a  +  Cj)2  +  (or  +  C^2  *  (1  -  v)  (a  *  C^cJ 

+  v{q*  Cx)3  Cj  «■  (a  x  Cj)4 
hg  -  +2k  («  x  Cp3  S1 

02 


(11  cont'd) 


h.y  ^  -k  (a  +  Cx)4 

h8  =  n2k[n2  -  4Sj2  -  (3  -*)<“  +  CpcJ*  <<*  +  C^2  Cj2 
+  2 i/(a  +  C.)3  C.  +  (a+  C  )4 
hg  =  k[-2n2  (a  +  Cj)  Sj  -  (a  +  Cj)  Sj3  -  2  (o  +  Cj)2  CjSj 

+  +  Cj)3  S1| 

h10  =  k|-2n2  (a  +  Ct)2  -  (a  *  Cj)2  Sj*  -  (1  +  r)  (a  +  Cj)3  Cj  | 
hn  =  -2k  (a*  C1)3S1 

h12  =  +k  <“  +  C!>4 
h  -  4<3  -  »>  “Si 

t2  -  4*1  +  “)“<“+  Cj) 
f3  -  ° 

t4  =  -p(a+  CjlCj  -j(l  -K)n2  -  Sj2 

*5  =  -<“+  C1)S. 
f6  -  (a  +  Ct)2 

*7  *  0 
f8  ’  '  °S1 

f9  =  (o  *  Cj)2  +  »<a  *  C,)  Cj 


G3 


g,  -n2(o+ Cj) -|(1  -  1-)  (a+ C^Sj2 +  |(1  -  k)  (q  +  Cj)2  Cj  (llcont'd) 

g2  =  "(I  -  ,)(o+  c1)2s1 

g3  =  +J-(1  -»)(<»+  Cj)3 

g4  =  +|-(3  -  »)  n  (or  +  Cj)  Sj 
g5  =  -f(l  +  k)  n  (a  +  Cj)2 
*6  “  *7  =  0 

gg  =  -n  (a  +  Cj)  Cj  -  vn  (a  +  Cj)2 

*9  =  g10  gll  g12  =  0 

Cj  -  (1  -  k2)  n[  -|-(a  +  Cj)2  CjM«*  Cj)3] 

c2  =  C3  =  0 

c4  =  -|-(1  -  v2){a+  Cj)3  Sj 

e5  =  (1  -  *'2)[-|-(2o+  Cj)  (a  +  Cj)3  +  la *  Cj)4] 

c6  =  °7  =  0 

Cg  =  (1  -  *-2> ["i-n2  +  C,)2>  (a*  Cj)3  Cj  *  (a  +  C,)4  | 

c9  =  -±-<l  -,2)(a+Cj)3Sj 

C10  =  +2"(1  '  '*><*«♦  C1><“*  C/ 

C11  =  C12  “  0 
*1  =  ‘2  “  a3  =  0 

a4  =  (1  -  M2)[y(2a  +  Cx)  (a  +  Cp  -  (o  +  C^2  +  yn2  J 
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0 


(11  concl'd) 


a9  =  <1  -  p2)|  -|-(2a+  0^(0+  Cj)  +  (a+  Cf] 

a10  =  *11  a12  =  0 

bl  =  U-*'2)[{(0’+C1)C12-(or+C1)2C1-i(2Qr+C1)(o+C1)C1-ioS12| 
b2  =  |-(1  -  a2)  |(or  +  C/  Sj  -  (20  +  Cj)  (o  +  Ct)  Sx  ) 
b3  =  -|il  -  A  (20  +  Cj)  (o  +  Cj)2 

b4  =  T(1  "  •’2>n<“*  ci>Si 
b5  ■  b6  ’  b7  '  0 

bg  =  (1  -  a2)  j+i-n  (o  +  Cj)  Cj  -  n  (o  +  C^2  | 
bS  ”  b10  =  bll  =  b12  =  0 


The  stability  equations  for  a  toroidal  shell  |Eqs.  ( 10)  |  consist  of  three 
linear  homogeneous  ordinary  differential  equations  with  variable  coefficients. 
The  unknowns  in  these  equations  are  the  three  displacement  components  uQty)  , 
v  {4>)  ,  and  w^)  .  For  a  complete  toroidal  shell,  the  boundary  conditions  are 
simply  conditions  of  periodicity  on  the  displacement  components;  hence,  «n(^)  , 
vn($)  ,  and  wr(^)  may  be  represented  by  Fourier  series  in  the  meridional 
coordinate  $>  .  Next,  in  order  to  pave  the  way  for  a  Fourier  series  analysis, 
the  coefficients  in  the  stability  equations  are  expressed  as  linear  combinations 
of  trigonome?  ic  functions.  Then,  with  the  aid  of  some  identities  which  will  be 
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derived,  it  is  possible  to  write  each  of  the  three  stability  equations  in  the 
following  form: 

OO  00 

1  |  |m  cos  <m«  ♦  2 

m=0  m=l 

The  braced  egressions  in  Eqs.  (12)  represent  homogeneous  linear  combinations 
of  the  Fourier  coefficients  used  in  the  expansions  for  un  »  vn  »  and  .  The 
trigonometric  functions  in  Eqs .  ( 12)  are  linearly  independent  for  0  £  ^  £  2  rr  ; 
hence,  each  of  the  braced  expressions  must  vanish.  Thus  the  problem  of 
solving  a  system  of  ordinary  differential  equation!  with  variable  coefficients  is 
reduced  to  that  of  solving  an  infinite  system  of  linear  homogeneous  algebraic 
equations.  A  matrix  iteration  technique  is  used  to  get  the  lowest  eigenvalue  of 
a  finite  system  of  equations  which  is  obtained  through  truncation  of  tho  infinite 
system  of  equations.  The  size  of  the  finite  system  of  equations  is  then  succes¬ 
sively  increased  until  no  significant  change  occurs  in  the  computed  eigenvalue. 
2.2  The  Solution 

In  order  to  facilitate  the  subsequent  Fourier  series  analysis,  we  now 
express  the  coefficients  in  the  stability  equations  as  linear  combinations  of  trig¬ 
onometric  functions.  This  is  effected  by  employing  the  well  known  trigonometric 
identities 

CaSb  *  llVb  '  Sa-b  1 
Ca°b  •  ilCa+b  +  Ca-bl 
B.S  *  ilCa-b  *  Ca+bl 


sin  (m^)  =  0  . 
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from  which  it  follows  that 


<“  *  cl>  C1 

=  T+ “ci  +TC2 

(13a) 

<«  +  Cj)  C2 

=  ici  +  ac2  +ic3 

(13b) 

<«  +  Cj)  Sj2 

_  Of  .  1  ft  ®  n  i  p 

2  4  U1  "  2  L2  ■  4  C3 

(13c) 

(“  +  Cl>  Sl3 

■  T^l  +  ls2  -  T“S3  -  Is 4 

(13d) 

(a  +  Cj)2  Cj 

-  “+(T+“2)ClTaC2+lC3 

(13e) 

<a  +  Cj)2  Cj2 

=  8"+X  +i<“Cl  +(t  +  T  )C2  +TC3  +  'TC4 

(13f) 

2  9 

(a  +  Cx) 

-  i.+  SL-  +  — C  -  —  C  -  — C  -  -i-C 

8  2  2  C1  2  C2  2  C3  T°4 

(13g) 

9 

9 

(a  +  Cj)2  CjSj 

=  Tsi  +(t+t)S2  +TS3  +Ts4 

(13b) 

(a  +  Cx)3 

*  i“+  °3  +(f +  30?)c1  +f“c2  +T°3 

(131) 

(a+  Cj)3  Cj 

=  !+I“2+(t“+“3)c1  +(t-P)c2 

+  —  aiP  4.  Ir 

(13j) 

4aC3+  8C4 

(a  +  Cx)3  Sx 

=  (-|-  a+  a3^  +(j+-|a2)s2  +4aS3  +TS4 

(13k) 

(o+  Cx)4 

=  J-+  3a2  +  a4  +  (3a  +  4a3)  Cx  +  (-|-+  3a2) C? 

+  aC3  +  lC4 

(131) 

(a  + 

=  oS1+TS2 

(13m) 

,  „  ,  „  2  a  3«  a  _  i  _ 

<“  +  Cl>  C±  =T^_C1+TC2+TC3  (131) 

(Of  +  C.)  C.S-  =  ~t^S-  +  -r-S_  +  ■4~S 
1  11  41  22  43 
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(13o) 


(13p) 


tQ+  ci>2  =T+“2  +  a*!  +f°2 

<“+ci>2si  "  (t  +  °2)si  +“S2  +  T®: 


(13q) 


Insertion  of  Eqs.  (13)  into  Eqs.  (10)  and  (11)  yields  the  following  form  of  the 
stability  equations: 

[<h10  +  c10x>  +  <hll  +  clix>  C1  +  <h12  +  C12X>  C2  +  <hl3  +  C13X>  C3]  V*> 

"  [h21  S1  +  h22  S2  +  h23  S3 1  Un^  "  I  (h30  +  h31  C1  +  h32  C2 
+  h33  C3]anW  +  [  <*41  +  C41X)  S1  +  (h42  +  C42X)  S2  +  <h43  +  C43X)  S3 
+  (h44  +  C44X)  S4]V^  +[  <h50  +  C50X)  +  *h51  +  C51X)  C1 
+  (h52  +  C52X)  C2  +  (h53  *  C53X)  C3  +  (h54  +  C54X)  C4|Vn(^ 

"  [h61  S1  +  h62  S2  +  h63  S3  +  h64  Vn(^  “  h  +  **71  C1  +  **72  °2 

+  *73  C3  +  **74  C4]  +  [(h80  +  C80X)  +  (h81  +  C81X)  C1 

+  (h82  +  C82X)  C2  +  (h83  +  C83X)  C3  +  (h84  +  C84X)  C4  I  wn(*> 

‘  hi  +  C91X^  S1  +  ^92  +  C92X^  S2  +  ^h93  +  C93X*  S3 
+  (h94  +  C94X)S4lWn(^  "  [(h10,0  +  C10,0X)  +  h10, 1  +  C10,1X)C1 
+  (h10,2  +  C10,2X)  C2  +  (h10,3  +  C10,3X)  C?  +  (h10,4  +  C10,4X)C4lWn(^ 
+  [hll,lSl  +  hll,2S2  +  hll,3S3  +  hll,4S4lWn^)  +  [h12,0 

+  h12, 1  C1  +  h12,2  C2  +  h12,3  C3  +  h12,4  C4lWn(^  =  0  <14a> 

fll  S1  ]un«*>  ‘  [f20  +  f21  CllUn(^  +  |(f40  4  a40X)  +  (f41  +  a41X)  C1 
+  <f42  +  a42X)  C2|VnW)  +  [f51Sl  +  f52S2lV*> 

I  f60  +  f61  C1  4  f62  C2l  Vn(^)  +  |f 8 1  Sl|Wn(V)  '  i  '*90  +  a90X) 

+  ^91  4  a91X^l  4  ^92  4  a92^  ^2!  Wr»^  ^  0 


G8 


[<*10  +  b10X>  +  <*11  +  bnX)Cl  +  <*12  +  b12X>S  +  <*13  +  b13X>C3|\<« 

*  I <*21  +  b21X)Sl  +  <*22  +  b22X)S2  +  <*23  +  b23X>  Ssl  “n»> 

'  I  <*30  +  b30X)  +  <*31  +  b31X)  C1  +  <*32  +  b32X>  C2 
+  <*33  +  b33X>  Si  +  |  <*41  +  b41X)  S1  +  <*42  +  b42X>  Si  V«> 

+  1*50  +  *51  C1  +  *52  C2  I V#  +  {<*80  +  b80X>  +  <*81  +  b81X»  C1 
+  <*82  +  b82X>SlWnW)  =  0  (14c) 

where 

hio  =  (^T i!)k  +  (1+ll')+  ^  ““ 

bll  =(|("T+T,’)+“2  +  "2]k  +  '4<1  +  *')  +  (1  +  3i>)a2Jn 

h12  =  |(“2  *  1f)k  +  (*  +Tl')“n 

b13  =  |*4(1  -■')*  +  7(1  +  ►)]  n 

b21  =(-i-+2o?)"k 

h22  =  +lank 

h23  =  ~4"nk 

! 

1 

j  h30  ~  -ank 

f 

b31  =  -{f*”*)14 

| 

j  h32  =  '0nk 

j  b33  =  -f1* 

'J 

! 
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41 


42 


*43 


44 


50 


'j«  ->  -b’  *  (> 


va 


a 


(15  cont’d) 


=  j  (1  -  v)  k 


-(*♦*■>]• 


k  .1(1  +  v)-±{l  +  3v}ot* 


-|-(1  -  *;)k  -|(1  +  V) 


51 


52 


53 


54 


61 


62 


63 


(t  +  \VY  (2  +  bV +  (i+  "V] k  +  t(1  +  v) 

+  (3  +  -|v  +  0?  jo? 

|^(11  +  9»0  +  (1  +  V)(?  +  2n2Jk  +  (3  +  -|*»») 

+  (4  +  1/)  c? 

||-(i  +  ^> +  (x  +lvK  +  in2|k  +T(1  +  ">  +  (3  +1*'K 
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-!■•*  +  |-<i  +  v) 


{t*  2q2)  Qk 

Hi+30?)k 


+|“k 


4* 


-(t-  3<^  +  “*)k 
-(3  +  4c?)  ok 


0 

Nl 
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Ns  “  ■Qk 


(15  cont'd) 


•>74  '  -Tk 

hso  “It*"7  +  "> +  n2|n2k  +  i<1  +  ^  +(t+  3,,  +  °2)“2 

hSl  =  I  *<3  ‘  ">  1,2,1  +  ^-t1  +  •'>  +  (■*  +  2f)  a2  |  <» 
hg2  «  -s-U  +  i/)  n\  +  (1  ♦  *)  *(y  *  3  v)a2 

h83  =l<1  +  ,,>“ 
h84  ■  h1  +  »> 

hji  -  l-i(-7  +  3i/)  +  va2  -  2n2|ak 
hg2  =  |<-3  +*>  +  (-1  +|-)o2-  n2  k 

1*93  *  -!*1  -,’>0,k 
*94  =  *5(1  '  ,,)k 

Vo  ■  h  -!■')-  (2 + b*  2i?y  -n2]k 

hio,  i  =  I  "T(n  +  9*4  -  (1  +  v)q2  -  **\ 

h10,2  ‘  4<1  +  ‘')-(1+f‘’)0?-n2  k 

;  "io.3  ■ -b +  *•** 

| 

! 

h10t  4  =  "?*,k 

|  *11.1  =(-!+2a2)ok 

!  hn.2  ■  ■§'**)* 

l 

I 

!  71 


hll,3 

hll,4 

h12,0 

h12,l 

h12,2 

h12,3 

h12,4 

C10 

C11 

C12 

C13 

C41 

C42 

C43 

c44 

C50 

C51 

c52 


(|-  +  3a2  +  a4)k 
(3  +  4a2)  ak 

(i+  3“2)k 

ak 

lk 

+(1  -  »/2)  (1  +  a2)  an 
+<1  -  *2)(f  ♦  fa2)n 
+(1  -  v2)an 

♦g-t1  -  */2)n 

-|  (!  -  v2)(-|  +  a2)  a 

-*<»-*(M«2) 

-fa  -  t?)o 

1  2, 

-16  ('  -  ”  > 

<>-*(**M 

(i  •  i'2,(i + \aiy 

<*  •  *'2)  (t  +  f“2) 


(15  cont’d) 
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(15  cont'd) 


„  -3/i  2. 

c53  8  ^  "  v  )° 


=  +_Lu  -  *2> 


'54  16 


'80 


'81 


'82 


'83 


'84 


'91 


'92 


'93 


(1  - 

A 

I* 

(*♦ 

2  1  2\  2 
o 

(1  - 

A\ 

'21 
i  4 

f  5o ? 

2  \ 

-  n  la 

(1  - 

A( 

1  + 

9  2 
2° 

1  2\ 

T  ) 

2, 

-  v  )a 

+T<> 

-  V 

■2> 

-(1  - 

A 

(i* 

1  2^ 
Ia) 

1  a 

_(1  - 

A 

Ik* 

3  2\ 
4«J 

I 

-1(1 
8  ' 

-  V 

2)a 

1  2 
4n 


„  -  JL/i  2. 

'94  “16  (1  “  v  > 


C10,0 

=  (i  -  A 

>(A 

+  la?  + 

4 

C10,l 

=  (1  -  v2) 

■ft 

*K)' 

C10,2 

-  (i  -  A 

F) 

C10,3 

-  *fr  - 

A<* 

°10,4 

“  tI5<1  ' 

■  ^2) 

fll 

’  +I«3  - 

p)  n 

f20  =  +2(1  +  v)m 
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(15  cont'd) 


f21  =  +  F<1  +  *'>“ 

f40  =  4(‘1+,’>+<1  --^l 
f41  *  * 
f42  =  +?(1 

f51  =  "  a 

f  =  _-L 
52  2 

t.  ■  Q  +  J) 
f61  =  +2° 
f62  ”  *2 
f81  = 

f9o  ■  +hi+v'*J 

f91  =  +(2  +  v)a 

f92  ■  +F<1  +  '> 

a40  "  (i  -  ^(-3-+i“2) 

a41  =  ‘I*1  *  ‘'2>“ 

8 42  ’  *7(1  -  "2> 

a90  -  Td  -  a2) 
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(15  cont'd) 


a9i =  +r(1.' 1,2)0 

a92  *  +T(1  -  *,2> 

*10 "  [i*1  -  >,)-"2|“ 

*n  =  <!  *  >'>(7*2)  • 1,2 

*12  -  Id  -»)» 

*13  =  4<X  ‘  1,1 

*30  =  (1  ‘  ‘’Jff  +  T)  0 

*31  ■ (1  ‘  ‘,)(i  +  |°2) 

*32  -  I'1 
g33  ~ 

g41  =  j(3  ~  v )  nor 

g4 2  =  T(3  “  V)  n 

*80  =  'I?*1  +  ">  +  1'“2|n 


g81 

g82 

g21 

g22 

g23 


-(1  +  2  m)  an 

l  /t 

-gU  +  v)  n 

-<■-<*  i) 

a 

-d  -  »)  y 
-|d  -  -) 
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(15  cont’d) 


*50  =  -<!  +  »)*(?  +  \) 
*51  =  *(1  +  x)  no 

*52  =  'T<1  +  *')  » 
bio  =  ‘T*1  *  2)a 

bu  =  -<1  -  x2)(|  +  202) 

b12  =  "f<1  "  2>a 

bi3 -  -h1  -  -2> 

b21  *  -J<2  -  *V 

b22  =  ‘T*1  -  ■'2»“ 
b23  *  0 

b41  ”  -|-(1  -  x2)  no 
b42  "  ~T{1  "  v  >n 

b8o =  ~(1 "  + °2)n 

b81  =  "l(1  "  l/2)an 
b82  =  -?<1  -  ^>n 
b30  =  “  *^)  U  +  o?)  a 

b3i  -  -<i--2)(l+!«2) 
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(15  concl'd) 


b32  =  -d  -  -V 

b33  ■  T*1  -  ■'2> 


For  each  value  of  n  ,  the  displacement  omponents  u^(  ip)  ,  vn(^)  »  and 
w n(^)  may  be  represented  by  the  Fourier  series: 


u„<« 

vn<« 

w  (ip) 
n'  ' 


^  Um  cos  m  ip  +  Um  sin 


m=0 


m=l 


00 


>  V  sin  +  >  V  cos  m ^ 

/ ,  m  m 


m=l 


m-~0 


Y  W  cos  m$  +  \  W  sin  mip 
/ .  m  m 


m=0 


m=l 


(16a) 


(16b) 


(16c) 


The  series  with  the  Fourier  coefficients  U  ,  V  ,  W  represent  a  buckling 

m  m  m 

mode  which  is  symmetric  about  the  plane  $  =  0,  *  (plane  A-A  in  Fig.  1), 

w 

whereas  the  series  with  the  Fourier  coefficients  U  ,  V  ,  W  represent  a 

m  m  m 

buckling  mode  which  is  antimetric  about  the  plane  *p  =0,  * .  After  insertion  of 
Eqs.  (16)  into  Eqs.  (14),  it  can  be  shown  that  the  resulting  stability  equations 
may  be  put  in  the  form: 


>  A  cos  m^  +  >  A  sin  m^ 
m  jL,m 


=  0 


(17a) 


m=0 


m=l 


y  II  sin  mip  +  >  B  cos  m$ 

ID  ^  111 


=  0 


(17b) 


m-1 


m-0 


y  C  cos  m*l>  +  \  C  sin  m*P 


-  0 


(17c) 


m=0 


m=l 
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where 


A  =  A  <U,  V  W  )  , 
m  m  m  m  m 


B  =  B  (U  ,  V  ,  W  )  , 
m  m '  m*  m*  m 


C  =  C  (U  ,  V  ,  W  )  , 
m  m '  m’  m’  m' 


(18a) 

(18b) 

(18d) 


B 

m 


=  A  (U  ,  V  , 
m  m  m 


=  B  (U  ,  V  , 
m  '  m  m 


=  C  (U  ,  V  , 
m  '  m  m 


W  )  (18d) 

m 

Wm)  (I8e) 
Wm)  .  (18f) 


Since  the  functions  sin  (m^)  and  cos  (m^)  in  Eqs.  (17)  are  linearly  independent 
for  0  <  ^  £  2  ir  ,  we  conclude  that 


A  =  B  =  C 

m  m  m 

»v 

A  =  B  =  C 

m  m  m 


0 

0 


(19a) 

(19b) 


From  Eqs.  (18  and  19)  we  see  that  the  Fourier  coefficients  U  ,  V  , 

m  m  m 

may  be  determined  from  a  set  of  equations  which  do  not  contain  U  ,  V  , 

mm 

W  .  This  means  that  a  toroidal  shell  under  uniform  external  pressure  can 
m 

buckle  into  a  mode  which  is  either  symmetric  or  antimetric  about  the  plane 
4>  =  0,  tt  ,  and  thus  these  modes  can  be  investigated  separately.  For  convenience, 
the  buckling  mode  which  is  symmetric  about  the  plane  ^  =  0,  v  is  called  Mode  A 
and  the  buckling  mode  which  is  antimetric  about  the  plane  =  0,  *  is  called 
Mode  B.  These  two  modes  are  considered  in  the  next  two  subsections. 

2.2.1  Mode  A 

For  the  buckling  mode  which  is  symmetric  about  the  plane  tp  =  0,  tt  ,  we 
let 
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(20a) 


00 

00 

u  (V)  = 

y  u 

cos 

mV 

=  y  u 

n 

Z,  m 

i—4  m 

m=0 

m=0 

00 

00 

ii 

C 

> 

y y 

Z*  m 

sin 

mV 

=  y  v 

Z*  m 

m=l 

m=l 

oo 

oo 

** 

II 

y w 

Z*  m 

cos 

i 

mV 

-  I  w» 

m=0 

m=0 

m 


(20b) 


(20c) 


where,  for  brevity,  we  have  used  the  notations  given  by  Eqs.  (3).  Substitution 
of  the  Fourier  series  expansions  for  the  displacement  components  |Eqs.  (20)J 
into  the  stability  equations  for  a  toroidal  shell  |Eqs.  (14)  J  yields 

[<h)0  +  C10X>  +  <hll  +  CaX>  C1  *  <h12  *  C12X)  C2 

00 

+  <h13  * UmCm-(h21Sl  +  h22S2 

m=0 

00 

+  h23  Ssl  S  m  Um  Sm  '  [h30  +  h31  C1  +  h32  C2 
m=l 

oo 

+  h33  I  Um  Cm  +  [  (h41  +  C41X>  S1  +  <h42  *  e42X>  S2 
m=0 

oo 

+  (h43  +  c43x>  S3  *  <h44  +  C44X>  S4]  I  Vm  Sm 

m=l 

+  Ko  v  C50X)  +  <h51  *  °5IX)  C1  +  <h52  +  C52X)  C2 

OO 

+  <h53  +  °53x>  C3  +  (h54  *  CS4X>  <=4]  2 .  m  Vm  Cm 

m~0 

(cont'd) 
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-  I  h61  S1  +  h62  S2  +  h63  S3  +  h64  S4|  I  “*  Vm  Sm 

m-1 


‘  [^O  +  C1  +  h72  C2  +  **73  C3  +  **74  C4l  J  ^  Vm  Ci 

m=0 

+  [<h80  +  °80X)  +  <h81  +  °81X)C1  +  (l>82  +  °82X)C2  +  <h83  +  C83X,C3 


+  <h84  +  C84X)  C4|  2  wm  cm  "  l(hCl  +  C91X)  S1  +  <h92  +  C92X)  Si 


m=0 


+  (h93  +  °93*>  S3  +  <h94  +  C94X>  S4!  I  m  Wm  Sm 

m=l 

•  l(Vo  +  Cl«,«x>  +  <hl«,l  +  C10,  Ix>  C1  +  <h10,2  +  C10,2x>  C2 

oo 

+  <h10.3  +  cl«.3x>  C3  *  <h10,4  +  C10,4X>  CJ  I  ”2  wm  Cm 

m=0 

OO 

+  |hll,lSl  +  hll,2S2+h11.3S3  +  hn,4S4l  1  m3wmSm 

m=l 


+  *  h12,0  +  h12,l  C1  +  h12,2  C2  +  h12,3  °3 


hio  aCaI  y  m4W  C  =  0 
12,4  4)  La  mm 

m=0 


(21a) 
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fll  S1  2  Um  Cm  ■  [f20  +  f21  Cl]  £  m  Um  Srr  +  [(f40  +  a40X) 
m=0 


m-1 


oo 


(f41  +  a4iX>  C1  +  (f42  ♦  a42X)  c2j  J  VmSm 


m=l 


[f51  S1  +  f52  S2]  2  m  Vm  C 
m=C 


-  [f60  +  f61  C1  +  <62  Si  2  m2  Vro  Sm 


m=l 


OO 


+  f81  S1  ^  Wm  Cm  [(f90  +  *90^  +  (f91  +  *91^  C1 


m=0 


+  <<92  T  a92X>  S]  I  mWmSm  '  0 


(21b) 


m=l 


[<*10  +  b10x>  +  <*11  +  bllX>  C1  +  <*12  +  b12X>  C2  +  >*13  +  b13X>  sl  2  Um  Cm 


m=0 


[(g21  +  b21X)  S1  +  (g22  +  b22X)  S2  +  g23  S3  ]  X  m  Um  Si 

m-1 

[*g30  +  b30X)  +  ^g31  +  b31X*  C1  +  ^g32  +  b32X*  C2 

OO 

+  <*33  +  b33X>  S)  2  ™2  Um  Cm  +  [<*41  +  b41X>  S1 


m=0 

OO 


+  <*42  +  b42X>  h\  1  Vm  Sm  +  [*50  +  *51  C1  +  *52  C2]  2  ■  Vm  Cm 


m=0 


m=0 


+  [<*80  +  b80X>  +  <*81  +  b81X>  C1  +  <*82  +  b82X>  S[  2  Wm  S,  =  0 


m^O 


(21c) 
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The  terms  in  Eqs.  (21)  consist  of  products  of  trigonometric  functions  and 


Fourier  series.  These  terms  can  be  put  into  the  form 

<0  00 
II  Im  Cm  •  or-  II  Im  Sm 

m=0  m=l 

by  use  of  some  transformation  equations  which  we  will  now  derive.  As  an 
example,  let  us  consider  a  term  of  the  form 

oo 

P  =  sr  I  m3  wm  sm  .  <£r<4)  .  (22) 

m=l 

We  use  the  trigonometric  identity 

2S  S.  =  C  .  -  C  . 
a  d  a-b  a+b 

to  write  P  as  the  sum 


2P  -  A-B 

where 


A  =  Y  m3  W  C 

m  m-r 

m=l 

00 

B  -  I  m3  Wm  Cm+r  • 
m=l 


(22a) 


(22b) 


(22c) 


In  Eq.  (22b),  we  let  q  =  m-r  .  Then 


A  =  Y  (q  +  r)3  W  C  , 
Zrf  '  q+r  q  ’ 

q--r+l 


or, 

oo  -1 

A  ■  I  <«  +  r>3  V  Cq  +  2  «>  +  r>3  Wq+r  Cq 

q=0  q=-r+l 


(22d) 


We  let  q  =  -p  in  the  second  ♦cnn  on  the  right-hand  side  of  Eq.  (22d).  Thus 


or, 


-1 


2  <1  +  r>3  V  cq  -  I  <r  -  P)3  wr-p  cp  . 


q=-r+l 


p=r-l 


-1 


r-1 


2  lfl*r|3VCq  *  ^(r-q)3W_.C. 


q=-r+l 


r-q  q 


(22e) 


q=l 


The  upper  limit  for  q  on  the  right-hand  side  of  Eq.  (22e)  can  be  changed  to  *> 
if  we  make  the  agreement  that  Fourier  coefficients  with  negative  subscripts 

vanish.  Hence 


-1 


y  (q  +  r)3  W  C 
Z,  Q+r  q 

q=-r+l 


I 

q~i 


w 


r-q 


C 

q 


(22f) 
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Insertion  of  Eq.  (22f)  into  Eq.  (22d)  results  in  the  following  expression  for  A 


w 

“  ^  (m  +  r) 


3wm+rCm  +  2  (r  -  m)3W„  „  C 


r-m  m 


m=0 


m=l 


+  r3W  C  -  r3W  C 
r  o  r  o 


<22g) 


3  3 

where  the  terms  tr  W  C  and  -r  W  C  were  added  to  the  right-hand  side 

r  o  r  o 


of  Eq.  (22d). 
But 


OO  OO 

y  (r  -  m)3 W  C  +  r3  W  C  =  S  (r  -  m)3W  C 
Z*  '  r-m  m  r  o  Z»  '  '  r-m  m 

m=l  m=0 

and 


-  r3  W  C  =  -  S  6  (m  +  r)3W  A  C 
r  o  Z*  mo'  '  m+r  m 

m=0 

where  6  is  the  Kronecker  delta: 
mo 

1 ,  m  =  0 
0  ,  m  *  0 

Hence ,  Eq.  (22g)  may  be  written  as 

oo 

A  ■  1  l*1  -  V”  +  r>3wm+r  +  <r  -  m>VJCm  '  <22h> 

m=0 
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By  proceeding  in  the  same  way,  we  can  write  B  in  Eq.  (22c)  as: 


B  =  V  <q  -  r)3  Wq_r  CQ  .  (22i) 

q=0 

We  now  insert  the  expressions  for  A  and  B  given  by  Eqs.  (22h  and  i)  into 
Eq.  (22a)  to  get 


00 

2P  =  y  |-(m  -  r)3  W  +  (1  -  6  )  (m  +  r)3W 

/L  1  m-r  '  mo'  '  '  m+r 

m=0 

Mr  -m)3Wr.m]Cm  •  <22j> 

Then,  with  the  notation 


+1,  m<r 
0,  m=r 
-1,  m>r 


we  can  write 


c 

mr 


-(m  -  r)3  W  „  +  (r  -  m)3  W  =  e  Im  -  rl3  W,  .  (22k) 

'  m-r  '  '  r-m  mr  |  m  -  rl  . 


Then  from  Eqs.  (22j  and  22k)  we  obtain,  finally, 


w 


S  >  mW  S 
m  m 

m=l 


"V 

-il  I 


€  lin 
mr 


rl3  W 


Im-rl 


m=0 


+  (1  -  6 _ )  (m  +  r)3  W  L  |c 

mo'  m+r  1  m 


(23) 
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Thus,  the  term  S  V  m**  W  S  has  been  expressed  in  the  form 
r  mm 

m=l 

00 

zL  I  |m  ~m  *  proceeding  in  the  same  way,  identities  may  be  obtained 
m=0 

which  enable  us  to  write  each  term  in  Eqs.  (21)  in  the  form 


I 

m=0, 1 


L  Cm  •  or'  I  I  L  S, 

m=l,  2 


These  identities  are: 


oo 

■  \ 

OO 

■  >  1 

2  C  >  W  C 
r  /_  mm 

I ! 

(1+6  )  W,  .  +  (1  -  6  )W  C 

mr  Im-rl  '  mo7  m+r  m 

(24a) 

m=0 

m-0 

OO 
■  % 

OO 

2S  >  V  S  = 
r  mm 

II 

+  €  V,  .  +  (1  -6  )  V 

mr  Im-rl  '  mo7  m+r 

k 

(24b) 

m=l 

m=0 

oo 

00 

1 

2S  V  W  C 

r  Za  mm 

II 

(1+6  )  W.  ,  -  W  S 

mr7  !m-rl  m+rJ  m 

(24c) 

m=0 

m^l 

OO 

■  H 

00 

2  C  >  V  S  = 
r  Z-r  mm 

I 

-e  V.  ,  +  V  Is 

mr  Im-rl  m+r |  m 

(24d) 

m=l 

m=i 

00 
■  s 

00 

2  C  >  m  V  C 
r  Zw  mm 

II 

im-ri  Vlm-rl  +  (1  -  6mo)<m  + 

r)  V  \ 

7  m+r  J 

C 

m 

nr=0 

m=0 

<24e) 

2  S  >  mW  S  = 
r  L->  mm 

00 

II 

tcmr  lm  -  rl  WIm-rl 

m-1 

m=0 

+  (1  -  6  )  (m  +  r)  W  |c 

mo7  '  7  m+r  J  i 


<24f) 


86 


OO  OO 

2Sr  I  mVmCm  =  I  Ilm  'rl 


Vlm-H  Vm+r  Sm  <24«> 


m-0 


m=l 


OO  00 

2  C  VmW  S  =  S  I- €  Im  -  rl  W  ,  +  (m  +  r)  W  S 
r  mm  \  mr  Im-rl  m+rJ  m 


m=l 


m^l 


(24h) 


2  C  V  m2W  C  =  >  Ilm  -  rl“  W.  . 

r  t__,  mm  \  Im-rl 

m-0  m-0 


Si 


+  (1  -  6  )  (m  +  r)2  W  |c 

'  mo  '  m+r  J  m 


(24i) 


2  S  >  m  V  S 
“  mm 

m-1 


OO 

-  SI 


*W  lm  -  rl  V Im-rl 


m^O 


+  (1  -  6  )  (m  +  r)2  V  |c 

'  mo'  '  '  m+r  J  m 


(24j) 


OO  °o 

2  Cr  I  Vm  Sm  *  S  Hmr  lm  *  1,2  Vlm-rl  +  <m  +  Vm+r  lS 


m-1 


m-1 


m 

(24k) 


2  sr  1  m2  Wm  Cm  ‘  1  ! Im  ‘  r|2  Wlm-rl  ‘  <m  +  r>'  Wm*rl  Sm  <241> 


m=0 


m=l 


ao 


2  C  >  mw  V  Cm 
in  m 

m=0 


■  I  I lm  -  r|3  Vlm-rl  +  «l  -  ‘J  <“  +  r>3  Vm+r  lC 


m=0 


m 
(24m) 


2  S  Ym3W  S  -  f\ 
r  ^  mm  Lu  \ 

m=l  m-0 


+  e  Im  -  rr  W.  . 
mr  Im-rl 


+  (1  -  <5  )  (m  +  r)3  vvlc 

mo  m+r  J  m 


(24n) 
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2Cr  m3W 


S 

m  m 


2  C 


m=l 

m=l 

y  m3V  C 

Z*  mm 

00 

-  II 

m=0 

m-1 

y  m4V  S 

Zj  mm 

ii 

M8 

m=l 

m=l 

Y  m4  W  C 
/,  mm 

oO 

■  21 

m=0 

m=0 

-  I  [  -£mr  -  r|3  Wlm-rl 


4  <"  +  r>3  Wm+r  |Sm 


(24o) 


cmrlm  -  rl  Vlm 


-ri  +  (m  +  r)4  Vm+r|S 


m 

(24q) 


‘W  C 
m+rl  m 


(24r) 


where 


mr 


mo 


mr 


1,  m=r 
0,  m^r  ,  r 

1,  m=0 

0,  m/0  , 

1,  m<r 
0,  m=r 
■1,  m>r  ,  r 


=  1  4 

At  •  •  f  t 


=  1,  ...  4 


(25) 

(26) 
(27) 


We  note  that  identities  involving  Um  may  be  obtained  by  replacing  in 
Eqs.  (24)  by  Um  . 

We  now  apply  the  transformations  given  by  Eqs.  (24)  to  the  stability  equa¬ 
tions  | Eqs.  (2 1) |.  As  a  result,  we  find  that  all  terms  in  each  of  Eqs.  (21)  may 
be  written  as  an  infinite  summation  over  the  same  trigonometric  function  and 
over  the  same  range  for  m  .  Hence,  in  each  of  Eqs.  (21),  we  may  collect  all 
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terms  under  the  same  summation  sign  and  arrive  at  the  following  form  of  the 
stability  equations  for  buckling  Mode  A: 


m=0 


2  2  <h10  ♦  C10X)  Um  ♦  J  (hlr  ♦  clrX)|(l  ♦  6mr)  U|m.r| 


r=l 


(1  -  !mo)0m.rl  '  2  h2r[£mr  lm  ’  rl  U|m-. 


r=l 


+  (1  -  6  )  (m  +  r)  U  I  -  2m2  hon  U 
'  mo7  '  7  m+rj  30  m 

3 

-  V  h_  f  lm  -  rl2  U,  ,  +  (1  -  6  )  (m  +  r)2  U  } 

3rl  Im-rl  '  mo' '  7  m+rJ 


r=l 


2  <h4r  +  V>|£mr  V-rl  +  <»  '  5JVr  1 


r=l 


2m  <h50  +  C50X)  Vm  +  2  <h5r  +  V>  [lm  ’  rl  V|m-rl 


r=l 


(1-6  J  (m  +  r)  V  I  -  Y  hc  [  €  lm  -  rl2  V, 

mo7  7  m+r  1  6r  l  mr  Im-i 


r=l 


Ml-«mo)(u>  +  r)2Vm+r|-2m3h,»Vm 
4 

-  2  h7r[lm  -  r|3  Vlm-rl  *  +  r>3  Vm+r  ] 


r=l 


+  2  <h80  +  V>  Wm  +  2  (h8r  +  VM'1  »  6mr>Wlm-ri 

+  *  6mo>  Wm+r  I 


r=l 
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-  X  (h9r  +  VX>  |cmr  '  rl  Wlm-rt  +  (1  ‘  6mo»  <m  +  r»  Wm+r  I 


r~l 


* 

-  2m2<hio,o  + ‘lo.o^m  '  I<hio,r  +  cio/>ilm-r' 


w 


Im-rl 


r=l 
4 


(1  -  «mo)(m  +  r)2Wm+r]+  £  hu,r(<mr  lm  -  rl3  W|m.r| 


r=l 


(i  -6mo)(m  +  r)3Wm+r|+  2m4h12>0Wm  +  j  hl2_r[im  -  ri4W|m_r| 


r=l 


+  d  -  «mo>  <»  +  r>4  wm+r] 


cos  (m  4)  -  0 


(28a) 


m=l 


1  I  'ill*1  +  «ml>  Ulm-ll  *  Vll*  2mf20Um  '  f2lhml lm  ‘  11  Ulm-ll 

+  <m  +  l>  Um+1 1  +  2  <f40  +  V>  Vm 

2 

+  S  ({.  +  a„  X)  l-f  V,  .  +  V  ! 

Z,  '  4r  4r  ;  I  mr  Im-rl  m+r  j 

r-1 
2 


2  *»!'“  -  rl  Vlm-rl  -  <“  +  r>VmJ-  2”2  f60  V 


m 


r~l 

2 


-  V  L  — €  lm  -  rl2  V,  -  +  (m  +  r)2  V 

6r  mr  Im-rl  '  '  m+r  J 


r~X 


f8 1 1  * 1  +  *ml>  Wlm-ll  -  Wm+ll"  2m  <f90  +  a90X>  Wm 


9r  +  a9r^  *mr  *m  r*  W|m-rl 


-2« 

r-l 

*<”*')  Wm+rl 


sin  (mV)  =  0 


(28b) 


00 


m=0 


S  2(g1A  +  b1AX)  U  +  y  (g,  +  b,  X)  (1  +  6  )U,  .  +  (1  -  6  )U 

Z*  I  610  10  '  m  Z-  lr  lr  '  mr'  Im-ri  '  mo'  m+rj 


r=l 


*  y  (g„  +  bn  X)  [€  Im  -  rl  U. 

Z  '®2r  2r  |  mr  Im-i 


r=l 


+  <l  -  6mo>  <m  *  r>  Um  J  ‘  2m2  (g30  *  b30X)  Um 
3 

-  I  <g3r  +  b3rX>  |lm  ‘  r|2  Ulm-rl  +  (1  ’  6mo>  <m  +  r»2  Um+r  I 


r=l 

2 


+  V  lgA  +  b„  X)  [  e  V.  .  +  (1  -6  )  V  I 

Z,  '&4r  4r  '  l  mr  Im-rl  mo'  m+r  I 

r=l 

2 

+  2m  g50  Vm  +  2  g5rl lm  '  rl  Vlm-rl  +  (1  '  U  (m  +  r>  Vm+r  I 


r=l 


2«80  +  b80X)Wm  *  I  <g8r  *  b8rX>l<‘  +6mr>Wlm-rl 


+  (1  -  6  )  W  A  | 

'  mo7  m+rJ 


r=l 

cos  (m^)  =  0 


(28c) 


In  order  for  Eqs.  (28)  to  be  identically  satisfied  for  all  values  of  ,  each 
of  the  braced  expressions  in  Eqs.  (28)  must  vanish.  This  yields 


z(m)u  +  y  z<m>u  +  y  z(m)u  +  Z(m)v  +y  z(m)v 

10  m  Z  2r  Llm-rl  Z  2r  m+r  Z  40  m  Z  Z5r  Vlm-rl 


r=l 


r=l 


r=  1 


y  z(m)v  -  z(n,)w  +  y  z(m)w 

Z  6r  m+r  z70  m  Z  z8r  Im-ri 


r=l 

4 

2 


r=l 


,(m) 


+  >  *Jr  wm+r  =  0  ’  m  =  «M.2 . ) 

r=l 


(29a) 
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<m)  (m)  (m)  (m) 

x10  Um  ^1  Ulm-ll  ^1  m+1  x40  Vm 


I 

r=l 


(m) 


+  >  v'm“'V  +  Y  v<m>V  +  v'“‘7  W 
Z  x5r  v Im-rl  Z  x6r  vm+r  x70  wm 


+  vlm) 


r=l 

2 


+  Y  x(m>  W  +  Y  x<®'  vJ 
Z  8r  Wlm-r!  Z  9r  Wi 


m+r 


0  ,  m  =  (1,2,...) 


r=l 

3 


r=l 

3 


(29b) 


r=l 


y<m)  u  +Y  y(m)  U  +  Y  y(m)  U  +  y(m)  V  +  Y  v(m)  V 
y10  m  Z  J2r  u  Im-rl  Z  y3r  m+r  y40  Vm  Z  y5r  Vlm-rl 

r=l 

2  2 

+  Y  v(m>  V  +  y(m)  W  +  Y  y(m)  W 
Z  y6r  vm+r  y70  Wm  Z  y8r  W Im-rl 


r=l 

2 


r=l 


2  y9r)wm+r  =  0  •  m  * 


r=l 


where 


(29c) 


rSm)  = 
ij 


i<ra> 

1J 


-(m) 
-  coz.. 

ij 


x; 


(m)  _ 


~(m) 
x..  ' 
V 


-(m) 

cox.. 

1J 


(30) 


(m) 

yij 


y(m) 

yij 


-(m) 

-  coy..  ' 


and 


•  -  i-.  1 


X  pa/Eh 


(31) 
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The  coefficients  on  the  right-hand  sides  of  Eqs.  (30)  are  given  by 


=  -1 

Z10  1 

2h10  -  2m2  h3ol 

(32) 

7<m)  -  - 
z2r  " 

*  amr>  hlr  '  lm  ’ 

rl  e  h 

mr 

*2r  ‘ 

lm 

-rl2 

h3r],  (r  =  1,2,3) 

S<m>  =  -a  -  6  ) 

3r  1  mo' 

hlr  -  (m  + 

r>  h2r  - 

(m  + 

r)2 

h3r  ] 

,  (r  =  1,2,3) 

Nl 

of 

II 

1 

2m  h50  - 

2m3  h?0  ] 

i<”>  =  - 

\e  h. 

i-  lm  -  rl  h_ 

-  lm  - 

rl2  e 

h* 

5r 

l  mr  4r 

5r 

mr 

6r 

- 

lm  -  rl3 

n 

1,2, 3, 4) 

-ZW  =  _n  -  5  ) 

z6r  '  mo' 

h4r  f  + 

r>  h5r  - 

(m  + 

D2 

h6r 

z(m) 

70 


-<m) 

z8r 


r<m) 

z9r 


-  (m  +  r)3  h?r]  ,  (r  =  1,2, 3, 4) 

~[2h80"  2m  h10,0  +  2m  h12, o] 

'll1  +6mr,h8r  -  lm-rl  emr  h9r  *  lm  '  r|2  h10,r 
+  !m  -  rl"'  €mr  hji  r  +  lm  -  r|4  h12> r]  ,  (r  =  1,2, 3, 4) 

•(1  ‘  6mo>  [h8r  -  <m  *  r>  h9r  '  <m  *  r>2  h10,r 


(m  +  r)Uhll  r  +  (m  +  r)  h12>rj  ,  (r  =  1,2, 3, 4) 


z(m)  =  2  c 
10  10 


z2r 


?(m) 

Z3r 


<1  +  *mr,clr’ 


(1-6  )  c,  , 

'  mo  lr 


(r  =  1,2,3) 
(r  =  1,  2,  3) 
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(32  cont'd) 


z(m)  =  2m  c 
Z40  C50 


5(m) 

Z5r 


z6r 


?(m) 

z70 


€mr  c4r  +  lm  '  rl  c5r  <r  =  1*2*3*4> 
<1-6mo>|c4r  +  <m  +  r>c5rl  <r  =  1-2-3-4: 


2  c80  ~  2m  C10, 0 


:<m)  . 
8r 


(1  +dmr>c8r  '  lm  '  rl  emr  c9r  *  lm  ‘  1 


;<m) 

x10 

-(m) 

X21 

-<m) 

*31 


-(m) 
x'  ' 
5r 


-<m) 

*70 

n(m) 

8r 


(1  -  Smo>|c8r  -  (m  +  r>c9r  *  <m  +  r>2  c10,, 

- 1"  2m  f  go  ] 

■|(1  +  aml)  £11  +  lm  ■  1 1  e  ml  f21 1 

-Hu  -  +  D  *21  i 

*(2  f40  -  2m2  f60  I 

'Hmrf4r  +  lm  -rl£5r+  lm  '  r|2  'mr  f6rl 
-(*4r  ’  <m  *  r>f5r  -  <m  +  r>2  f6r  I  *  <r  ’  1 
-|-2m  f90  I 

-[(1  +  amr)arlf8r  +  lm  ’  rl  f  mr  f9r  I  •  <* 

'Kl  f8r  -  <m  +  r>  f9r  I  >  <r  =  ».  2) 
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10, r  (r  =  1*2*3*4) 

] » (r  =  1,2,3, 4) 


,  (r  =  1,2) 

,2) 

•  =  1,  2) 


(32  cont'd) 


$(m)  _ 
X31 

0 

*(m)  _ 

x40 

2  a40 

*<.”>  = 
*>r 

"fmr  &4r 

C(m)  _ 
x6r 

a4r  <r  =  1»2) 

ii 

-2m  a9Q 

0<m>  _ 
x8r 

,m  -  rl  emr  a9r 

A(m) 
*9r  " 

-(m  +  r)  agr  (r  = 

v(m)  = 
y10 

-|2  *10  -  2m2  *30 

— (m) 

V  = 

>+6mr>*lr- 

^<m>  = 
y3r 

-<!  -amo»l*lr  - 

-(m) 
y40  = 

-|2m  *50 1 

-<m) 

y^r  = 

~|€mr  *4r  +  lm  * 

— (m) 
y6r  = 

*<*  -Smo»[*4r  + 

— (m) 
y70  “ 

"  2|®8ol 

-0*0 
nr'  - 

-  <»  +  V  *8r  (r 

.2  1 

-  Im  -  rl  g3r]  , 


*5r 


y9r 


“  -  ®mo>*8r  (r  “  1'2> 


'(m)  _ 
^10 


2  bio  -  2m  b30 


<i  =  1,2,3) 

1,2,3) 
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-<m)  s 
y2r 

=  (1  +  6  ) 
'  mr' 

blr 

-  Im 

"  rl  fmr  b2r  " 

Im  -  rl2 

b3r  »  <r  =  1.2,3) 

(32  concl'd) 

0<m>  = 
y3r 

=  (1  -6  ) 
'  mo' 

lblr 

-  (m 

+  r)^  -  (m 

+  r>V 

)  (r  =  1,2,3) 

<<> 

ii 

=  0 

«(m) 
y5r  = 

=  e  b. 
mr  4r 

(r  = 

1,2) 

£<m)  . 
y6r 

b. 

4r 

(r  = 

1,2) 

ii 

?  o 

o. 

2  b80 

C  -  <1+6mr>b8r  <r=1'2> 

C  '  <l  -  5  mo»  b8r  <r  -  ™ 


The  coefficients  a  ,  b  ,  c  ,  f  ,  g  ,  and  h  in  Eqs.  (32)  are  given  by  Eqs.  (15). 

By  letting  m  take  on  the  values  m  =  0,1,2 _ in  Eqa.  (29),  we  obtain  an 

infinite  system  of  algebraic  equations  in  which  the  unknowns  are  the  Fourier 

coefficients  U  V  and  W  .  The  coefficients  in  this  system  of  equa- 
m  m  m 

tions  are  shown  in  Table  1.  Hence,  the  system  of  ordinary  differential  equations 

with  variable  coefficients  |Eqs.  (10  and  11)  J  has  been  replaced  by  an  infinite 

system  of  linear  homogeneous  algebraic  equations  |Eqs.  (29),  or,  see  Table  1 J 

for  U  ,  V  ,  and  W 
m  m  m 

Using  matrix  notation,  we  rewrite  Eqs.  (29)  as: 


E  V  -  wS  V  =  0 


(33) 


where  |  R  j  and  |  S  j  are  square  matrices  formed  by  the  coefficients 
(SjJm)  »  y|jm)  .  *{”*)  and  .  y}™*  ,  zj™*)  .  respectively;  j  V  |  is  a 
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column  vector  formed  by  the  unknown  Fourier  coefficients  U  ,  V  ,  W 

m  m  m 

The  elements  r^  ^  and  s^  ^  of  [R]  and  [S]  can  be  obtained  from  Table  1  and 
Eqs.  (30,  32,  and  15).  For  example,  the  diagonal  element  s_  in  the  5^  row 


of  [S]  may  be  expressed  in  terms  of  the  geometric  parameters  ~  and  -  ,  and 

a  a> 

the  number  of  circumferential  waves  n  ,  as  follows: 


*5,5 


Table  1  &  Eqs. 


(30)  +  z*1* 

S  z70  Z82 


Eqs.  (32)  -  2  hg0  -  hg2  +  hQO  +  2hin  n  +  h 


92 


-  h  _  o  h  _  h 

‘10,0  “10,2  11,2  12,0  n12, 2 


Eqs.  (15) 


_1 

12' 


-  3  -  v 


3  o  2 
-v-2o 


y 


where  the  reference  above  the  long  equal  sign  indicates  the  means  by  which  the 
right-hand  side  was  obtained.  Similarly  the  element  rg  g  =  z^q  +  . 

For  fixed  values  of  the  geometric  parameters  a  =  ^  and  k  =  »  and 

for  a  fixed  value  of  the  number  of  circumferential  waves  n ,  the  coefficients 
z<”>  ,  x<"»  ,  and  in  Eqs.  (33)  depend  only  on  the  load  parameter  X. 

The  objective  of  the  stability  analysis  is  to  determine  the  lowest  value  of  X  (or 
equivalently,  the  highest  value  of  w  =  l/X)  for  which  Eqs.  (33)  admit  a  non¬ 
trivial  solution.  This  value  of  X  is  called  an  eigenvalue  of  the  system  of  equa¬ 
tions  and  the  associated  nontrivial  solution  is  called  an  eigenvector.  The 

components  of  the  eigenvector  are  the  Fourier  coefficients  U  ,  V  ,  and  W 

mm  m 

These  components  are  only  determined  up  to  a  scalar  factor.  The  displacement 
functions 


OO 

u  (ip)  =  >  U  cos 
n  Z*  m 

m=0 


m  ip 


(34a) 


97 


Coefficients  of:  U  W 

u. 

V, 

w, 

u. 

v0 

w„ 

u. 

V„ 

W, 

u 

o  o 

1 

1 

i 

2 

2 

2 

3 

3 

3 

r  8 

9 

10 

11 

12 

13 

i4  n 

,<0) 

10 

y(°> 

>70 

y<0) 
y  21 

y(°; 

y51 

(0) 

y81 

y(0) 

>22 

y<0) 

>52 

y(0) 

>62 

y<°) 

>23 

! 

,<0) 

TO 

z<°) 

Z70 

z(°) 

Z21 

zi°) 

Z51 

z<°) 

Z81 

■8 

2(0) 

52 

7<°> 

62 

z<°) 

23 

z<°> 

Z53 

z(°) 

83 

z(0) 

Z54 

(°) 

*84 

(1) 

21 

xd> 

81 

xd> 

10 

xd>  +  xd> 

40  52 

xd>  +  xd) 

70  82 

4? 

X(D 

61 

xd> 

X91 

xd> 

62 

xd) 

92 

(1) 

2’. 

yd> 

y81 

yd)+yd> 
y10  y22 

.,d)+yd) 

y40  >52 

yd)+yd) 

y  70  y  82 

vd>+vd> 

>23  >31 

(-> 

>61 

yd) 

91 

yd> 

>32 

yd) 

>62 

d) 

>92 

yd) 

>33 

(1) 

21 

zd> 

81 

zd>+zd) 
z10  z22 

zd)+zd> 

z40  Z52 

*d>  +  zd> 

z70  z82 

zd)  +  ld> 

Z23  31 

z(1> 

53 

+  z<') 
Z61 

z(1> 

83 

*4? 

z(1> 

32 

zd)+zd) 
Z54  Z62 

7U)+Zd) 

84  92 

zd) 

Z33 

zd> 

z63 

z(1) 

z93 

x<2) 

82 

x(2) 

21 

x<2) 

X51 

x<2) 

81 

x<2) 

10 

x(2> 

40 

\'2> 

70 

x(2) 

31 

x<2> 

X61 

X<2> 

91 

x(2) 

X62 

x<2> 

92 

(2) 

22 

y(2) 

y82 

y(2)+y(2) 
y21  y23 

y(2) 

y51 

y(2) 

y81 

y<2> 

y10 

y(2) 

>40 

y<2) 

>70 

y(2) 

y31 

v« 

>61 

y(2) 

y91 

y(2) 

>32 

y(2) 

>62 

y(2) 

>92 

\2) 

,<2) 

.(2)  12) 

,(2)  .  ,(2) 

.(2)  (2) 

,(2) 

,(2) 

j.  .(2) 

J2) 

.J2, 

,(2) 

J2) 

J2) 

.(2) 

.(2) 

.(2) 

.0 

H)  .(4) 


THIS  PATTERN  REPEATS 
FOR  m  -  5 


(m)  ,.<m) 


(m)  Cm) 


Jra) 

y(m) 

„(m) 

*21 

X51 

xbl 

L.(m) 

v.(m> 

i 


(m)  _(m) 


Tabic  1 

MODE  A  STABILITY  EQUATIONS 


00 


(34b) 


Vn  ^  =  2  Vm  sin 
m=l 

00 

{Ip)  =  y  W  cos  uf  (34c) 

X*  xu 

m=0 
.  i 

corresponding  to  the  eigenvalue  X  are  called  the  eigenfunctions  or  mode  shapes. 
In  order  to  effect  a  solution,  we  truncate  the  infinite  set  of  equations  |  Eqs.  (33)  J 
and  employ  a  matrix  iteration  technique  [Ref.  (26))  to  determine,  from  Rayleigh's 
quotient,  the  lowest  eigenvalue  X  of  the  resulting  finite  equation  system.  Next, 
the  size  of  the  equation  system  is  increased  and  the  lowest  eigenvalue  of  the  new 
system  is  determined.  This  procedure  is  repeated  until  the  successive  values 
of  X  have  stabilized  to  the  value  X  .  Then  by  varying  r. ,  we  obtain  a  set  of 

vl 

these  values  of  X  ,  one  for  each  n  ,  and  the  solution  to  the  buckling  problem 
cr 

is  given  by  the  minimum  value  of  X  in  this  set.  This  procedure  was  pro- 

vl 

grammed  for  the  IBM  7094  digital  computer  and  the  numerical  results  are 
presented  in  the  next  chapter.  At  this  point  we  may  remark  that,  although  a 
mathematical  proof  of  convergence  has  not  been  presented  here,  the  numerical 
results  indicate  convergence  of  the  eigenvalues  within  a  wide  range  of  values  of 
geometric  parameters. 

2.2.2  Mode  B 

For  the  buckling  mode  which  is  antimetric  about  the  plane  $  =  0,  * 

(plane  A-A  in  Fig.  1),  the  displacement  components  ufl  (if>)  ,  (^)  ,  and  wr(^) 

are  represented  by  the  Fourier  series: 

CO  oo 

un  <W  "  I  8in  <”«  =  2  Sm  (35a) 

m=l  m=l 
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oo 


(35b) 


v  (4>)  =  Y  V  cos  mit  =  y  V  C 

n  'r  ZL  m  ^  m  m 

m=0  m=0 

OO  OO 

wn  m  *  I  sln  =  I  Sm  •  <3Sc> 

m=l  m=l 

A  detailed  analysis  of  buckling  in  the  mode  which  is  symmetric  about  the  plane 
=  0,  ir  has  been  given  in  the  preceding  subsection.  Corresponding  equations 
for  the  antimetric  mode  can  be  obtained  in  the  same  way  and  therefore  a  detailed 
analysis  for  this  buckling  mode  will  not  be  presented  here. 

The  stability  equations  for  the  buckling  mode  which  is  antimetric  about  the 
plane  =  0,  *  are 
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where 


;(m) 

Zi3 

-  :(m)  _ ,  ,r(m) 

■  u£ij 

(37a) 

~(m) 

y 

-  ;(m)  _  ,  j(m) 

-  «a  -“{ij 

(37b) 

r 

-(m) 

y« 

_  ;(m)  .  Jrn) 

^ij  ^ij 

(37c) 

The  coefficients  on  the  right-hand  sides  of  Eqs.  (37)  are  given  by 

C  *  *|2  h10  -  2n>2  h3o]  <38> 

=  -|-fmr  hlr  +  lm  *  rl  h2r  +  fmr  lm  *  1,2  h3rl  «r  =  1>2-3> 

*3“*  *  -jhlr  -  (m  +  rjh^  -  (m  +  r)2h3r]  (r  =  1,2,3) 

*  -|-2”*>50+  2n,3h7ol 

"  'I'1  +  Smr>  h4r  *  fmr  ,m  *  rl  h5r  ‘  ‘  r|2  »6r 

€mr  lm  ‘  r|3  Nr  ]  (r  °  1>2»3»4> 

I  *  *£>  =  -(-h^  -  <m  +  r>  h5r  4  <n>  +  r>2  h6r  +  <>”  +  r)3  hr,)  (r  -  1.2,3, 4) 


4/4  4 


(38  cont'd) 


T(m) 

fc70 

Ttm) 

8r 


2  .  „4 


T(m) 

9r 


-2h80-2m‘h10>0  +  2m"h12,0 


-hmr  h8r  +  lffi  *  rl  Nr  +  f  mr  lm  '  r|2  h10,r  *  lm  *  r|3  hll,r 
*  fmrlm  *  r|4h12.rl  <r  =  1-2-3-4> 

*|h8r  -  <m  +  r)  h9r  -  <m  +  r)2  h10.r  +  (m  +  r)3  hll.r 

+  (m  +  r)4h12>rj  <r  =  1»2*3*4) 


£<m)  =  2c 


10 


e(m) 

i2r 


-Wclr  <r  =  1*2*3> 


C(m) 

3r 


=  clr  (r  =  1,2,3) 


-  -2m  c 


40 


50 


?(m) 

5r 


(1  +  t  )  c.  +  e  lm  -  rl  c_  (r  =  1,2, 3, 4) 
'  mr7  4r  mr  5r  '  ' 


-c4r  -  (m  +  r)c5r  (r  =  1,2,3, 4) 


C(m) 

70 


2  C80  ”  2m  C10,0 


t(m) 

8r 


-cmr  c8r  *  lm  ‘  rl  c9r  +  fmr  lm  ‘  rl  c10.r  <r  =  1-2-3>4> 


= 

9r 


c8r-<m  +  r)Co^  -  (m  +  r)  c1ftT.  (r  =  1,2, 3, 4) 


"9r 


lOr 


*£>  -  -12m  f 


‘20 


T<m)  = 
21 


*mi  fn +  ta  - 11  '2i 


3?'  -<1-Smo>|'ll  +  <"1+1>'21 
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(38  cont’d) 


F(m) 

40 


-I2  f40  -  2m  f6Q 


T5r'  '  -[<1+Smr>f4r-  'mr  -  rl  f6r  -  In,  -  rl2  I 


T6?’  =  -(1  -  «mo>lf4r  -  <m  +  r)f5r  '  <m  *  ^  f6r 


T(m)  =  -  2m  f 
*  70  21,1  *90 


C  *  -Kl  W  f8r  +  -  rl  f9r  1  <r  =  l-2> 

C*  =  -*1  *  4n,o>Klf8r  +  <m  +  r>  f9r  I  <r  *  ^ 


;(m) 

MO 


=  0 


?(m)  _ 


21 


=  0 


i<m>  =  o 

*31 


|(m) 

40 


=  2a 


40 


f(m) 

5r 


<1  +  6mr>a4r  <r  =  X’2> 


;(m)  _ 

'6r 


(r  =  1.2) 


£<*n) 

*70 


=  2m  a, 


90 


|(m) 

*8r 


Im  -  ri  a9r  (r  =  1,2) 


|(m) 

9r 


rj<m> 

10 


<*  ’  4mo>  (“  +  r>  V  <r  *  J'2» 


=  -[2  *10  ‘  20,2  *30 1 


Tj<m) 

2r 


"I"  fmr  «lr  +  lm  “  rl  *2r  +  €mr  ,m  *  rl  *3r 


6r  I  <r  *  *-2> 
(r  =  1,2) 


(r  =  1,2,3) 
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f<m>  = 
3r 


=  -|glr  -  <m  +  r)g2r  -  (m  +  r)2  g3rJ  (r  =  1,2,3) 


(38  concl'd) 


Tf<m>  = 

'40 


-  -k.m  g. 


rT(m>  = 
5r 


rj(m>  = 
6r 


=  'I*1  +  6n,r>  *4r  +  fmr  lm  ’  rl  g5r  I  <r  =  1’2) 
=  -|-g4r  -  (m  +  r)g5r|  <r  =  1,2) 


If*1”!  = 
'70 


=  -|2*U 


7f<m>  = 
8r 


~~€mrg8rj  1,Z^ 


nim)  = 

9r 


-  ‘hr  I  <r  *  ^ 


C’  '  2b10  -  2n,2b30 


2r 


-t _ b,„  Im  -  rl  b2r  +  fmr  in?  -  rl  bar  (r  =  1,2,3) 


mr  lr 


=  b.  -  (m  +  r)b9  -  (m  +  r)  b«  (r  =  1,2,3) 


=  <1  +  6mr)b4r  <r  *  ll2) 


«<“) 

6r 


-b4r  (r  =  1.2) 


"  €mr  b8r  <r  =  1'2) 


(j<m> 

9r 


bgr  <r  -  1,  2) 


The  coefficients  a  ,  b  ,  c  ,  f  ,  g  ,  and  h  In  Eqs.  (38)  are  given  by  Eqs.  (15). 
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By  letting  m  take  on  the  values  m  -  0,1,2 .  in  Eqs.  (36),  we  obtain 

an  infinite  system  of  linear  homogenous  algebraic  equations  in  which  the 

unknowns  are  the  Fourier  coefficients  IJ  ,  "V  ,  and  W  .  The  coefficients 

mm  m 

in  this  system  of  equations  are  shown  in  Table  2.  Using  matrix  notation,  we 
rewrite  Eqs.  (36)  as: 


5l|v|  -“[s]jv|  =  I  o} 


(39) 


where  J  R  j  and  |  S  j  are  square  matrices  formed  by  the  coefficients 

(I}®’  ,  5  <”>  .  £<“>)  and  (f<“>  .  if<“>  .  respectively;  |  v|is  a 

column  vector  formed  by  the  unknown  Fourier  coefficients  U  ,  V  ,  and  W 

J  mm  m 

The  buckling  load  is  obtained  by  the  method  discussed  in  the  preceding 
subsection  and  numerical  results  are  presented  in  the  next  chapter. 

3.  Stability  Equations  for  a  Sphere 

The  stability  equations  for  axially  symmetric  buckling  of  a  sphere  subject 
to  gas  pressure  are  oDtained  by  setting  u  =  n  -  a  -  0  in  the  stability  equa¬ 
tions  for  a  torus  jEqs.  (1C)|  .  As  a  result,  Eq.  (10c)  which  corresponds  to 
equilibrium  in  the  circumferential  direction  is  identically  satisfied.  And,  with 
the  aid  of  the  trigonometric  identities  given  in  Eqs.  (13). the  non-zero  coefficients 
in  Eqs.  (10a  and  10b)  are: 


h4  -’kli(S  - -)S2  *  i'l  -»>S4|-<1+»)|tS2  *  Js4]  -  C14  k[<2  -P)T1 

+  T3]  -  (J  <  V)  Tj  (40a) 

f  2!  I  (40b) 

kj(i  +  t')  +  Tx  |  +  (1  +  v) 
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where 


q  =  ~(1  -  p2)X  , 


(41) 


and 


Tj  =  tan  i 


(42) 
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With  a  trivial  change  in  the  independent  variable 


.  IT 

*  =  *  +  2  ’ 

we  arrive  at  the  following  form  of  the  stability  equations  for  symmetric  buckling 
of  a  sphere  under  hydrostatic  pressure: 

(1  +v  )[v  cot0  +  v  +  2wJ  +  k[-v  (2  -  v  +  eot20)  cot  0  +  v  (1  +  v  +  cot20) 

-  2  cot  0  v  -  v  +  w  (2  -  v  J-  cot2  0 )  cot  0  -  w  (1  +  v  +  cot2  0 )  (43a) 

+  2  cot  0  w  +  w  J  +  q  J  v  cot  0  +  v  +  4w  +  w  cot  0  +  w  J  =  0 

~(„  +  cot2 0 )  v  +  v  cot  0  +  v  +  w  (1  +  V)  -  q  (v  -  w)  =  0  .  (43b) 

These  equations  are  identical  to  Eqs.  (VU-76a)  and  (VII-76b)  of  Ref.  (13) 
provided  that,  as  was  done  here,  the  transverse  shear  force  is  omitted 
from  the  equation  of  equilibrium  in  the  meridional  diiection  [Eq.  (VII-76a) 
of  Ref.  (13)]. 


VII 


TOROIDAL  SHELL  UNDER  EXTERNAL  PRESSURE  - 
NUMERICAL  RESULTS 


Numerical  results  for  the  buckling  of  a  complete  toroidal  shell  under 
uniform  external  pressure  are  presented  in  this  chapter.  All  numerical 
results  are  for  a  value  of  Poisson's  ratio,  v  =  0.3  . 

1.  Numerical  Results 

The  procedure  used  to  obtain  the  lowest  eigenvalue  will  now  be  illus¬ 
trated  by  means  of  a  typical  example  Let  us  consider  a  toroidal  shell 
whose  geometric  parameters  are  a/h  =  100  and  b/a  =  4  ,  and  which 
buckles  in  Mode  A  with  n  =  2  circumferential  waves.  From  these  values, 
all  of  the  coefficients  r^  £  >  sk  £  *n  infinite  system  of  stability  equations 
[Eqs.  (VI  -  33)]  may  be  determined.  A  matrix  iteration  technique  is  then 
used  to  get  the  lowest  eigenvalue  X  of  a  finite  system  of  equations  which  is 
obtained  through  truncation  of  the  infinite  system  of  equations.  The  size  of 
the  finite  system  of  equations  is  determined  from  the  number  of  harmonics  M 
used  in  the  series  expansions  for  the  displacement  components  un($)  ,  vn($)  . 
and  w ^(ip) .  That  is,  in  Eqs.  (VI  -  20),  the  Fourier  index  m  takes  on  the 
values  m  =  0  ,  1  ,  2  , . . .  ,  M .  The  results  obtained  from  the  matrix  iteration 
method  are  shown  in  Table  lfor  a  system  of  equations  corresponding  to  M  =  14. 
From  Table  1 ,  we  see  that  four  place  accuracy  in  the  eigenvalue  X  =  pa/Eh 
was  achieved  after  eight  iterations.  Next,  by  assigning  a  sequence  of  values 
to  M  ,  we  successively  increase  the  size  of  the  system  of  equations  until  no 
significant  change  occurs  in  the  computed  eigenvalue  X .  The  resuits  so 
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obtained  are  shown  in  Table  2.  We  see  from  Table  2  that  the  sequence  of  eigen- 

-3  .  .  . 

values  X  converged  to  the  value  X  =  0.1746  *  10  .By  proceeding  in  this 
way,  we  can  determine  the  eigenvalues  X  corresponding  to  different  integer 
values  of  the  number  of  circumferential  waves  n .  Such  results  are  given  in 
Table  3  which  shows  that  the  lowest  eigenvalue  occurred  at  n  -  2  .  To  com¬ 
plete  the  analysis  of  the  toroidal  shell  with  a/h  -  100  and  b/a  -  4  ,  it  is 
necessary  to  consider  also  buckling  of  the  shell  in  Mode  B  and  in  the  axially 
symmetric  mode  (n  =  0) .  The  eigenvalues  for  these  modes  can  be  obtained 
by  the  same  procedure  used  for  Mode  A. 


Table  1 

RESULTS  OF  MATRIX  ITERATION 


Mode  A:  a/h  =  100,  b/a  =4,  n  =  2,  M 

=  14 

Iteration 

Iteration 

1 

7.6551 

9 

0. 1746 

2 

1.0454 

10 

0.  1746 

3 

0. 1973 

11 

0. 1746 

4 

0.  1763 

12 

0. 1746 

5 

0.  1752 

1  o 
io 

0. 1746 

6 

0.  1747 

14 

0. 1746 

7 

0.  1746 

15 

0. 1746 

8 

0.  1746 
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Table  2 


EFFECT  OF  SIZE  OF  MATRIX  ON  COMPUTED 
BUCKLING  PRESSURE 


Mode  A:  a/h  =  100,  b/a  -  4, 

n  =  2 

M 

1h*103 

M 

ft*1"3 

1 

5.  0164 

13 

0.  1746 

2 

1.  1389 

14 

0. 1746 

3 

0.4549 

15 

0. 1746 

4 

0.2762 

16 

0.  1746 

5 

0.2136 

17 

0.  1746 

6 

0.  1840 

18 

0.  1746 

7 

0.  1773 

19 

0.  1746 

8 

0.  1749 

20 

0.  1746 

9 

0.  1746 

21 

0.  1746 

10 

0.  1746 

22 

0.  1746 

11 

0. 1746 

23 

C.  1746 

12 

0.  1746 

24 

0. 1746 

Table  3 

VARIATION  OF  BUCKLING  PRESSURE  WITH 
NUMBER  OF  CIRCUMFERENTIAL  WAVES 


Mode  A:  a/h  = 

100,  b/a  = 

=  4 

n 

E2-*  103 

Eh 

n 

ft  * 103 

1 

0. 6978 

6 

0.  2878 

2 

0.  1746 

8 

0.3873 

3 

0.  1923 

10 

0.5167 

4 

0.  2179 

12 

0.  6714 

5 

0. 2494 

18 

1. 2427  j 
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The  effect  of  the  size  of  the  system  of  equations  on  the  computed  eigen¬ 
value  is  illustrated  in  Table  4  for  other  values  of  the  geometric  parameters 
a/h  and  b/a.  Inspection  of  Table  4  reveals  that  the  number  of  harmonics  M 
required  for  the  same  degree  of  accuracy  in  the  eigenvalue  increases  with 
increasing  a/h  and  decreases  with  increasing  b/a. 

For  all  cases  considered,  the  eigenvalues  for  the  axially  symmetric  buck¬ 
ling  mode  were  higher  than  the  eigenvalues  for  the  asymmetric  buckling  modes. 
For  the  same  value  of  n ,  the  eigenvalues  corresponding  to  the  two  asymmetric 
buckling  modes  (i.  e.  ,  Modes  A  and  B)  were  always  close  to  each  other,  and,  as 
can  be  seen  from  Table  5,  the  eigenvalues  for  Mode  A  were  sometimes  higher 
and  sometimes  lower  than  the  eigenvalues  for  Mode  B.  In  all  cases  investi¬ 
gated,  the  lowest  eigenvalue  occurred  at  n  =  2  for  both  asymmetric  buckling 
modes.  Some  mode  shapes  for  n  =  2  are  given  in  Figs.  1  and  2. 

For  the  limiting  case  of  a  sphere  (b/a  —  0)  under  external  pressure,  the 
classical  solution  p  =  1.21  Eh  /a  was  here  reproduced  numerically  for  both 
the  asymmetric  and  axially  symmetric  buckling  modes.  Results  for  the  limit¬ 
ing  case  of  axially  symmetric  buckling  (b/a  —  °°)  can  be  compared  to  the 
critical  load  for  an  infinitely  long  cylinder  under  external  pressure.  However, 
in  order  to  make  the  analyses  comparable,  we  first  have  to  modify  the  cylinder 
analysis  through  deletion  of  the  transverse  shear  force  from  the  in-surface 
equilibrium  equation.  Once  this  has  been  done,  the  results  from  the  two 
analyses  are  identical. 
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Table  4 


EFFECT  OF  SIZE  OF  MATRIX  ON  COMPUTED  BUCKLING  PRESSURE 


Mode  A,  n  =  2 

a/h=500,  b/a=2 

a/h=100,  b/a=2 

a/h=500,  b/a=8 

a/h=T00,  b/a=8 

M 

It x  lo4 

ft**3 

o 

r-H 

X 

«j|.c 

ft*”3 

6 

3.9804 

0.5176 

0.  3104 

0. 1157 

7 

2.3602 

0.3886 

0.  2195 

0. 1134 

8 

1.5499 

0.3203 

0.  1626 

0.  1132 

9 

1.0051 

0. 2949 

0.  1403 

0. 1132 

10 

0.7757 

0.2848 

0.  1283 

0. 1132 

11 

0.5824 

0. 2824 

0.  1239 

0.  1132 

12 

0.4828 

0.2813 

0.  1222 

0. 1132 

13 

0.4072 

0.2812 

0.  1217 

0.  1132 

14 

0. 3652 

0.2811 

0.  1216 

0.  1132 

15 

0. 3367 

0.  2811 

0.  1216 

0.  1132 

16 

0. 3201 

0.2811 

0.  1216 

0. 1132 

17 

0.3108 

0.  2811 

0.  1216 

0. 1132 

18 

0. 3052 

0.  2811 

0.  1216 

0. 1132 

19 

0.3029 

0. 2811 

0.  1216 

0. 1132 

20 

0.3015 

0.2811 

0.  1216 

0. 1132 

21 

0. 3010 

0.  2811 

0.  1216 

0. 1132 

22 

0.3008 

0. 2811 

0.  1216 

0. 1132 

23 

0.3007 

0.2811 

0.  1216 

0. 1132 

24 

0.3007 

0.  2811 

0.  1216 

0. 1132 

30 

0. 3007 

1 _ 
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Table  5 


COMPARISON  Or  RESULTS  FOR 
MODE  A  AND  MODE  B 


pa 
Eh  X 

104 

Mode  A 

Mode  B 

0.441 

0.440 

0.301 

0.301 

0.191 

0.  192 

0.  121 

0.  122 

0.068 

0.069 
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Fig.  1  Mode  Shapes,  a/h  -  100,  n  =  2  ,  Mode  A 
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Fig.  2  Mode  Shapes,  a/h  =  100 ,  n  -  2  ,  Mode  B 


1. 1  Buckling  Curves 

The  accuracy  of  the  numerical  results,  of  course,  depends  on  the  use  of 
a  sufficient  number  of  terms  in  the  Fourier  series  expansions  for  the  displace¬ 
ment  components.  On  the  other  hand,  the  computer  time  increases  rapidly 
with  an  increasing  number  of  terms.  The  convergence  of  the  method,  there¬ 
fore,  was  explored  through  calculation  of  buckling  loads  for  fixed  shell  param¬ 
eters  and  a  successively  increasing  number  of  terms.  By  use  of  these 
exploratory  calculations  it  was  possible  to  establish,  as  a  function  of  the 
geometrical  parameters,  the  number  of  terms  needed  for  1%,  or  better, 
accuracy  in  the  final  results. 

The  computed  critical  values  of  the  external  pressures  are  shown  in 
Fig.  3.  From  Fig.  3  we  see,  as  expected,  that  the  critical  pressure  p  de¬ 
creases  with  increasing  a/h  and  increases  with  increasing  a/b . 

1.2  Rigid  Body  Modes 

Let  us  denote  the  coefficients  of  the  unknowns  U0 ,  V  ,  Uj ,  ,  and  Wj 

in  the  i**1  row  of  [S]  or  [SJ  by  a.,b.  ,c.  ,d.,  and  e^  respectively. 
Then,  from  Eqs.  (VI  -  15,  29,  30,  and  32)  or  Eqs.  (VI  -  15,  36  through  38). 
we  can  obtain,  for  each  value  of  i ,  the  following  relations: 

(i)  Mode  A,  n  =  0 


(ii)  Mode  A,  ~  -  1 


(iii)  Mode  B,  n  =  0 


a.  - - c. 

i  a  i 


ai  =  ei  "  di 


di  =  -  ei 
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Fig.  3  Buckling  Coefficients  for  Toroidal  Shells  Under  Uniform  External 
Pressure 
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(iv)  Mode  B,  n  =  1 


bj  =  -a(d.  +  e.)  -  c.  .  (4) 

With  these  relations,  we  can  show  that  the  equations  which  govern  the  stability 
of  a  toroidal  shell  (Los.  (VI  -  29)  for  Mode  A  or  Eqs.  (VI  -  36)  for  Mode  B] 
admit  the  following  nontrivial  solutions: 

(i)  Mode  A,  n  =  0 

X  =  0  (5a) 

u  =  r  =  a(a  +  cos  ip)  (5b) 

v  =  0  (5c) 

w  =  0  (5d) 


(ii)  Mode  A,  n  =  1 

X  =  0 
u  -  sin  9 
v  =  sin  ip  cos  9 
w  =  -  cos  ip  cos  9 


(6a) 

(6b) 

(6c) 

(6d) 


(iii)  Mode  B,  n  =  0 

X  =  0 
u  =  0 

V  =  COS  Ip 
w  =  sin  ip 

(iv)  Mode  B,  n  =  1 

X  =  0 

u  =  sin  ip  pm  9 
v  =  (1  +  a  cos  ip)  cos  9 
w  =  or  sin  ip  cos  9  . 


(7a) 

(7b) 

(7c) 

(7d) 

(8a) 

(8b) 

(8c) 

(8d) 
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Nov  these  solutions  which  occur  at  a  zero  value  of  the  eigenvalue  X,  are 
recognized  as  being  the  following  rigid  body  modes  for  a  toroidal  shell: 

Case  (i)  Rotation  about  the  axis  of  revolution  (Mode  A,  n  =  0) 

Case  (ii)  Translation  in  a  plane  normal  to  the  axis  of  revolution 
(Mode  A,  n  =  1) . 

Case  (iii)  Translation  along  the  axis  of  revolution  (Mode  B,  n  =  0) 

Case  (iv)  Rotation  about  an  axis  normal  to  the  axis  of  revolution 
(Mode  B,  n  -  1) 

[Two  more  rigid  body  motions  are  obtained  from  Eqs.  (6  and  8)  through  re¬ 
placement  of  9  by  9  +  ir/z  1 

Since  the  proposal  here  is  to  obtain  buckling  loads,  it  is  necessary  in  each 
of  these  cases  to  determine  the  lowest  non-zero  value  of  X  for  which  the 
stability  equations  admit  a  nontrivial  solution  and,  of  course,  this  value  of  X 
was  used  to  arrive  at  the  numerical  results  already  presented. 

2.  Comparison  With  Test  Results 

A  comparison  of  the  results  of  the  present  theory  for  complete  toroidal 
shells  and  available  results  from  tests  conducted  at  Lockheed  Missiles  &  Space 
Company  is  given  in  Table  6a.  From  Table  6a,  we  see  that  test  and  theory 
agreed  to  within  10%. 

We  note  that  the  present  results  for  complete  toroidal  shells  can  be  used 
to  predict  buckling  pressures  for  partial  toroidal  shells  which  are  simply 
supported  on  the  meridional  edges  9  =  0  ,  9  =  p  .  For  example,  the  solution 
for  a  simply  supported  shell  with  fi  =  n  ,  2  tt/3  ,  tt/2  ,  . . .  can  be  obtained 
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COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


from  the  solution  of  a  complete  toroidal  shell  wMtih-buckles  in  n  =  2, 3, 4, . . . 
circumferential  waves.  A  180°  toroidal  shell  was  tested  at  LMSC.  A  com¬ 
parison  between  test  and  theory  is  shown  in  Table  6b.  Again,  the  agreement 
was  within  10%. 

The  infinitely  long  cylinder  (b/a  =  «)  and  the  sphere  (b/a  =  0)  repre¬ 
sent  limiting  cases  of  a  toroidal  shell.  For  external  pressure  loading,  it  is 
well  known  that  the  correlation  between  theory  and  experimental  results  is 
reasonably  good  for  the  cylinder  and  quite  poor  for  the  sphere.  Now  the  cor¬ 
relation  between  the  present  theoretical  results  and  the  few  available  test 
results  shown  in  Table  6a  was  quite  good.  However,  the  test  results  were 
for  a  slender  torus  (b/a  =  8)  ,  and  it  should  be  pointed  out  that  the  correla¬ 
tion  might  not  be  as  satisfactory  for  smaller  values  of  b/a  . 

3.  Comparison  with  Previous  Investigation 

The  only  previous  investigation  of  the  stability  of  a  toroidal  shell  under 
external  pressure  was  performed  by  Machnig  (Refs.  9  and  10).  In  the  first  of 
his  papers  (Ref.  9) ,  Machnig  studied  both  axially  symmetric  and  asj'mmetric 
buckling  modes  and  concluded,  contrary  to  the  present  results,  that  the  former 
buckling  mode  gives  the  smallest  critical  pressure.  In  his  more  recent  paper 
(Ref.  10),  Machnig  considers  only  the  axially  symmetric  mode.  A  comparison 
between  Machnig' s  results  and  the  results  of  the  present  analysis  is  given  in 
Fig.  4. 
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VIII 

FREE  VIBRATIONS  OF  PRESTRESSED  SHELLS  OF  REVOLUTION 

The  components  of  the  applied  loads  per  unit  area  of  the  shell’ s  middle  sur¬ 
face  are  denoted  by  p  A9,<t>)  .  pfl(0,  <£)  .  and  P  (0.0)  as  shown  in  Fig.  IV-1. 

The  shell  is  in  a  state  of  equilibrium ,  called  the  prestressed  state  of  equilibrium , 
under  the  action  of  the  applied  loads  and  the  resulting  membrane  stress  resultants 
N  (0,*)  ,  N0O(0>4>)  »  ’  311(1  ‘  The  effects  of  deformation 

of  the  prestressed  shell  will  be  neglected.  The  additional  quantities  that  develop 
as  the  shell  vibrates  about  its  prestressed  state  of  equilibrium  are  denoted  by 
(see  Figs,  m-2  and  3)  N^(0,<M)  ,  •  •  ,Q0(0,<M) ,  .  .  ,  M^0(0,<M) ,  •  •  ,  w(0,0,t) . 
These  additional  or  incremental  quantities  are  considered  to  be  infinitesimal. 

The  mass  per  unit  area  of  the  shell’s  middle  surface  is  denoted  by  p  . 

1.  Basic  Equations  for  Shells  of  Revolution 

The  equations  of  motion  for  a  differential  element  of  the  vibrating  shell  are 

2  2  2  2 

obtained  by  addition  of  the  inertia  terms  -  rr^pd  v/9t  ,  -rr^nd  u/9t  ,  and 
-rr^B  w/9t  to  the  left-hand  sides  of  the  equations  of  force  equilibrium 
derived  for  the  stability  problem  of  a  shell  of  revolution  (Eqs.  (IV  -  5)] : 

2fi  ■  <rN/  +  riN^  -  riN„ cos  *  +  f-%  -  rN«o"» '  rIWV 
+  ■•rl<N90“zl)’  -  <rVo“Zl>'  -  rlN9*o"z2  008 

9^v 

+  Sphrrlpzu,e  -  rV  ^2  ■  0  <la> 


124 


2F2  ‘  <rV  +  rlN'o  +  rlN(*  008  *  4  {-rlQ9  sin  *  '  rlNW“«  sin  * 

■  rlN?o"<f,  sin  *>  +  1  -rlK»o"*l  C0S  *  +  <rN»o“z2»'  4rl<NS*o"z2>’> 


,2 

+  6  .rr, p  cv ,  -  rr,u  ^  =  0 
ph  l*z  <£  r  2 


(lb) 


2F3  =  -riN9  sln  *  -  rtV  riQ'«  '  (rV  ‘  (rN«o‘V  '  ri(%o"/ 

-ri<W'  -  «r,W  +  |rN«eo"zi  -  riNWo“z2  8in  *1 


4  vrip*(‘>+c0)-rri,iil,2  =  0 

(7 1 

(1C) 

2mj 

■  <rM<./  4  rlMi  4  riM9«A  COS  1>  -  -Tl%  =  ° 

(Id) 

2M2 

■  -<rV  -  rlM«^  +  rlM9  COS  *  4  rrlQ*  ‘  ° 

vie) 

The  rotations  ,  o>zl ,  and  c*>z2  are  defined  in  Eqs.  (Ill  -  8)  through 

Eq.  (HI  -  10)  and  6^  is  defined  in  Eq.  (IV  -4). 

The  elastic  laws  which  relate  the  incremental  stress  resultants  N  ,  ,  .  .  , 

<t> 

Q. »  •  •  »  M  . .  to  the  incremental  displacement  components  are  given  by 
o  <pd 

Eqs.  (IE  -  18). 

2.  Shells  of  Revolution  Under  Axially  Symmetric  Loads 

From  Eqs.  (1),  we  see  that  the  vibrations  of  a  shell  of  revolution  are 
governed  by  a  system  of  partial  differential  equations  with  variable  coeffi¬ 
cients.  For  the  case  of  axially  symmetric  loading  the  coefficients  in  these 
equations  are  independent  of  the  circumferential  coordinate  0  and  the  time  t . 
Consequently,  it  is  possible  to  separate  the  space  and  time  variables  and  thus 
replace  the  system  of  partial  differential  equations  by  a  system  of  ordinary 
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differential  equations.  Such  a  separation  of  variables  is  effected  by  means  of 
the  following  Fourier  series  representation  for  the  incremental  displacement 


components: 


ito  t 

»,<M)  =  y  \(4>)  sin  (nfl)  e  n 


y(e,<t>,  t) 


oo 

-  2  V 


(<f>)  cos  (nd)  e 


OO  , 

XT'  10>nt 

w(0,tf>,t)  =  y  wn(0)  cos  (n0)  e 


Then  by  proceeding  in  the  same  way  as  in  Chapter  V ,  we  obtain  a  set  of  three 
ordinary  differential  equations  for  the  three  displacement  components  , 

vn(<f>),  and  w n(4>)  •  For  a  given  value  of  n,  these  equations  are 

<hl  +  °1>U„  4  <h2  +  «2>6n  +  <h3  +  C3>fin  +  <h4  4  °4)Vn  4  (h5  +  °5)i'n 
+  (h6  +  °6>  Vn  +  (h7  +C7>Vn+(h8  +  C8  +  68)Wn  +  <h9  4  °9>  *n 
+  <h10  +  °10>  #n  4  (hll  +  cll>Wn  4  <h12  +  °12>  %  =  0  (3< 


<fl  +  al>“n  +  <f2  4  a2,iln  +  <£3  +  a3>fin  4  (t4  +  a4  +  Vv„  4  <f5  4  a5)(rn 

4  lf6  *  4  (f7  +  a7)Vn  +  (f8  4  a8)wn  4  (f9  4  a9>*n  4  (fl«  4  a10>"n 

.«  .  •••  -  ,  •••* 

4<fll4all'wn4<f124a12>wn  =  0 
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tel  +  bj  +b1)%  +  te2  +  b2)un  +  (g3  +  b3)un  +  (g4  *  b4)vn  +  (g5  +  b5)vB 
+  (g6  +  b6>  \  +  (g7  +  b7)V„  +  <g8  +  b8>  Wn  +  teg  +  b9>  wn  +  <gw  +  bw)  w, 

+  teii  +  bii)wn  +  tei2  +  b12)»„  =  0 


(3c) 


where 


+  TT^lU) 


2 


(4) 


and  all  other  coefficients  in  Eqs.  (3)  are  given  by  Eqs.  (V  -  9).  For  convenience, 

a  subscript  n  was  deleted  from  a,..  ,h,  a^,  cgl  and  u>  in  Eqs.  (3) 

and  (4).  We  note  that  the  functions  z.(<p) ,  b{<p)  ,  and  c(<p)  in  Eqs.  (3)  depend 

on  the  prestress  quantities  N  .n  ,  N  A  ,  and  p  ,  whereas  the  coefficients  in 

<pU  ffU  z 

Eq.  (4)  depend  on  the  frequency  co  . 

3.  Free  Vibrations  of  a  Prestressed  Toroidal  Shell 

The  equations  governing  the  free  vibrations  of  a  prestressed  toroidal  shell 
are  obtained  through  specialization  of  the  equations  for  a  general  shell  of  revo¬ 
lution  [Eqs.  (3)].  We  consider  the  case  in  which  the  prestress  is  due  to  a 
uniform  pressure  p ;  p  is  positive  for  external  pressure  and  negative  for 
internal  pressure.  The  notation  for  a  toroidal  shell  with  a  circular  meridian 
is  shown  in  Fig.  VI-1. 

It  can  be  shown  that  a  toroidal  shell  under  an  initial  uniform  pressure 
can  vibrate  in  a  mode  which  is  either  symmetric  or  antimetric  about  the 
plane  $  =  0  ,  ir  (Plane  A-A  in  Fig.  VI-1).  The  symmetric  mode  is  called 
Mode  A,  and  the  antimetric  mode  is  called  Mode  B.  These  two  modes  are 
considered  separately  in  the  next  two  subsections. 
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3. 1  Mode  A 


For  the  vibration  mode  which  is  symmetric  about  the  plane  $  =  0  ,  it , 
we  let 


00  oo 

U(4>)  =  Y  U  cos  nvp  =  y.  U  C 
n  m  Y  m  m 


m  =  0 


m  =0 


(5a) 


oo  oo 

v(^)  =  V  V  sin  m$  =  £  V  S 

n  **  m  T  ,  ™ 

m  =  1  m  =  1 


m  m 


(5b) 


ay 

w  (*)  =  S  W  cos  mip  =  2  w  c 

n  m  „  m  m 


m  =  0 


m  =0 


(5c) 


where,  for  brevity,  we  have  used  the  notations  given  by  Eqs.  (VI  -  3).  Then 
by  inserting  Eqs.  (5)  into  Eqs.  (3)  and  proceeding  in  the  same  way  as  in 
Chapter  VI,  we  obtain  the  following  form  of  the  vibration  equations: 


z(m)U  + 
Z10  m 

1 

m-r  | 

♦ 1 4? 

r  =  1 

r  =  1 

i  4?v, 

m-r  | 

+  i  4f 

r  =  1 

r  =  1 

4?>Wm  + 

i  4”>w, 

m-r  | 

+  i  4? 

m+r 


m+r 


m+r 


0  ,  m  (0,1,2,...) 


r  =  1 


r  =  1 


(6a) 


x(ra)u  +  x^m^U  +  x^m^U  +  x^mV  +  V  x/inV  +  Y  x^mN 
x10  um  21  |m-l|  31  um+l  x40  m  Z  5r  |m-r|  Z  6r 


2 
i 

r  =  1 


2 

l 

r  =  1 


2  6 

+  x^mNv  +  Y  x^m^W  +  Y  x^m^W  =  0  m  =  fl  2  1 

x70  wra  Z  8r  |  m-r|  Z  9r  m+r  ’  (i,^,...) 

r  =  1  “  ' 


2 

l 

r=  1 


(6b) 


m+r 
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v<m>u  +  Y  y*m^U  +  Y*  y^m^U  +  y^m^V  +  Y  y^m^V 

‘10  m  Z  y2r  u|m-rj  2*  y3r  m+r  *40  Vm  2*  y5r  Im-rj 
r=l  1  1  -  1 


2 

l 

r=  1 


r  =  1 
2 


2 

I 

r  =  1 


+  V  v(m)V  +  v(m)W  +  Y  v(m)W  +  Y  v(m)W 
Z  y6r  m+r  y70  m  Z  y8r  w|  m-r|  Z  y9r  W 


r  =  l 


2 

l 

T-  1 


=  0 


where 


■r  ■  v  -  §  (c  -  -n 

4m)  ■  r-  ±i*r-*ri 

(m)  *(m)  _  1  [-(ro)  _  x*(m)l 

yij  yij  O  iyij  yij  I 


m+r 

m  =  (0,1,2,...)  (6c) 

(7a) 

(7b) 

(7c) 


and 


A  =  ^ 
A  Eh 


(8) 


and 


S2  = 


2  2 
^a  i o* 

Eh 


(?) 


The  coefficients  z|^  ,  ,  yj™^ ,  zj^  ,  ,  and  yjj*^  in  Eqs.  (7)  are 


given  by  Eqs.  (VI  -  32).  The  remaining  nonzero  coefficients  in  Eqs.  (7)  are 

J(m)  =  2c 
z70  c80 


*r  -  <»  +  «mr)E8r  •  (r  “  1 , 2 , 3 , 4) 


“  <*  -  Vsr  ■  <r  =  1 ' 2  •  3 ' 4> 


?m)  =  2a 
x40  40 


(10a) 

(10b) 

(10c) 

(lOd) 
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where 


;<m>  =  -c  a, 
or  mr  4r 


J(m) 

x6r 


(r  =  1,2) 

(r  =  1,2) 


y10 


J(m) 

y2r 


;<m) 

y3r 


"  (1  +  6mr>Blr 
"  (1  -  6mO>Blr 


(r  =  1,2,3) 


(r  =  1,2,3) 


580  =  "(1  "  ^(f  +  3a2  +  a4) 

cgl  =  -(1  -  w2)(3a  +  4a3) 

s82  =  -a-^l  +  3®2) 

S83  =  -(1-«'2)a 
®84  =  '  p  )  g 

sio =  (i  ■  -Mi a + “3) 

6n  =  (1  •  ^)(! +  3q2) 

B12  =  (1  -  VH\  a) 

big  =  (1  -  v2)  \ 


(lOe) 

(lOf) 

(10g) 

(lOh) 

(101) 

(lla) 

(lib) 

(lie) 

(lid) 

(lie) 

(Ilf) 

dig) 

(llb) 

(lit) 
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(llj) 


*40  *  <»  -  +  “2) 

*41  =  (1  -  v2)2 a 

hi  *  (» -  ■’*>  i 


(llk) 

(lll) 


In  Eqs.  (10)  and  (11),  we  have  used  the  notations  given  by  Eqs.  (VI  -  6),  (VI  -  25), 
(VI  -26),  and  (VI  -27). 

By  letting  m  take  on  the  values  m  =  0,1, 2,...,  in  Eqs.  KZ) ,  we  obtain 

an  infinite  system  of  algebraic  equations  in  which  the  unknowns  are  the  Fourier 

coefficients  U  ,  V  ,  and  W  .  The  coefficients  in  this  system  of  equations 
mm  m 

are  shown  in  Table  VI-1.  Using  matrix  notation,  we  rewrite  Eqs.  (6)  as 

[R)  {V}  -  (i)(S  -  XT)  {V}  -  {0}  (12) 

where  [R]  and  [S  -  XT]  are  square  matrices  formed  by  the  coefficients 


respectively;  {V}  is  a  column  vector  formed  by  the  unknown  Fourier  coef¬ 
ficients  U  ,  V  ,  and  W  .  The  elements  of  the  [R]  ,  [S]  ,  and  [T] 
mm  m 

matrices  can  be  obtained  from  Table  VI-1  and  Eqs.  (7),  (10),  (11),  and  (VI  -  15). 
3.2  Mode  B 

For  the  vibration  mode  which  is  antimetric  about  the  plane  f  =  0  ,w 
(plane  A-A  in  Fig.  VI-1),  the  displacement  components  are  represented  by  the 
Fourier  series: 

oo  OO 

u  ($)  =  U  sin  m^  =  X  US  (13a) 

n  _  _  i  m  ,  mm 

m  —  i  m  —  l 
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(13b) 


00 


v  (i£>)  =  V  cos  mill  = 
n  “  m 


00 


m =0 


n£o 


c 

m  m 


uo  00 

wM)  -  V  W  sin  mip  =  ^  W  S 

n'T/  _*■*  m  T  m  m 


m  =  1 


m  =  1 


(13c) 


The  vibration  equations  for  Mode  B  are 


3  3 

z^U  +  y  zlm)  U,  +  7  zim^  U 

10  m  Z*  2r  I  m-r  I  Z<  3r  m+r 

1  r  =  1 


r  =  l 


+  z<”>v  +  y  2<m> v,  ,  +  f  z'm)v  „ 

40  m  Zj  5r  m-r  6r  m+r 

r = 1  1  1  r=l 


+  z<”>w  +  y  z<m>w 

70  m  Zj  8r 


r  =  1 


+  V  2(m)W 
!  m-r  I  Z  9r  m+r 

1  '  r  =  1 


=  0  ,  m  =  (1,2,...)  (14a) 


2  i 

;(m)u  +x(m)u  +  x(m)U  +  x(m)v  +  Y  x(m)V  +  Y  x(m)v 
x10  m  21  U|m-l|  31  Um+1  40  m  Z  5r  lm-r|  Z  6r  m+r 

r=l  '  r=l 


+  +  y  x(m)W  +  y  x^m)W  =  0 ,  m  =  (0,1,2,...)  (14b) 

70  m  Z  8r  |  m-r  |  L*  9r  m+r 
r  =  1 


I  m_r  | 


r  =  l 


v(m)U  +  V  v(m)U  +  y  v(m)U  +  y(m)V  +  y  y(m)V 
y10  Um  Z  y2r  |  m-r  |  Z  y3r  >n+r  y40  Vm  Z  y5r 

r=  1 


|  m-r 


r  =  1 


2 

> 

r  =  1 


m-r 


2 

+  y  v^m^v  +  v^m^w  +  y  v^m^w  +  y  y^m^w 

Z  y6r  Vm+r  y70  m  Z  y8r  w|  m-r  I  Z  y9r  i 

r  =  l  r  =  1  r  =  1 

ra  =  (1,2,...) 


2 

l 

r  =  1 


I  m-r  I  ^  J9r  "m+r 
r  =  1 


=  0 


(14c) 
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where 
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L  ij 
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1J  1 
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- 
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f(m) 
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{(”) ,  i 

sij  '] 

(m)  p(m)  - 
Lj  ’  £ij  ’  4 

,  and  n!m)  in  Eqs. 
ij 

(15)  are 

-38). 

.  The  remaining  nonzero  coefficients  in  Eqs. 

(15)  are: 

J(m) 

41 70 

=  2S80 

(16a) 

?(m) 

&8r 

"€mrc8r 

> 

(r  =  1 

,2,3,4) 

(16b) 

J(m)  =  £ 

&9r  8r 


>>  =  25 
5  40  40 
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-  e  b. 
mr  lr 


,  (r  =  1 , 2 , 3 , 4) 
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"  <1  -  ‘mo^r  •  <r  -  1  •  2> 


*(m)  _  r 
^r  blr 


,  (r  =  1 , 2 , 3) 

,  (r  =  1 , 2 , 3) 


(16c) 
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Using  matrix  notation,  we  rewrite  Eqs.  (14)  as: 


[ft]  M  -(£)[S  -  XT){V}  =  {0} 


(17) 


where  [ft]  and  [§  -  XT]  are  square  matrices  formed  by  the  coefficients 
f*(m)  *(m)  i(m)l  ,  fr(m)  _  ^(m)  -(m)  _  -(m)  -(m)  _  -(m)l 

l*ij  ’^ij  *£ij  I  “  l«ij  *«ij  >T?ij  %  *£ij  ^ij  I* 
respectively;  {V}  is  a  column  vector  formed  by  the  unknown  Fourier  coef¬ 
ficients  U  ,  V  ,  and  W  .  The  elements  of  the  [R]  ,  [S]  ,  and  [T] 
mm  m 

matrices  can  be  obtained  from  Table  Vl-2  and  Eqs.  (11),  (15),  (16),  and 


(Vi  - 15). 

5.3  Eigenvalues  and  Eigenfunctions 

The  eigenvalues  and  corresponding  eigenfunctions  for  Modes  A  and  B  can 
be  obtained  from  Eq.  (12)  and  Eq.  (17),  respectively.  For  a  given  value  of 
the  pressure  parameter  X ,  the  lowest  eigenvalue  (2  for  which  these  equa¬ 
tions  admit  a  nontrivial  solution  may  be  obtained  by  the  same  method  as  was 
used  in  Chapter  VH  for  the  stability  analysis  of  a  toroidal  shell.  The  corre¬ 
sponding  eigenfunctions  or  mode  shapes  can  be  determined  after  substitution 
of  the  computed  eigenvalue  in  Eqs.  (12)  and  (17). 
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IX 


CONCLUDING  REMARKS 


This  work  presents  a  theoretical  investigation  of  the  stability  of  thin  shells 
of  revolution.  Stability  equations  are  derived  for  a  shell  of  revolution  under 
general  loading  conditions.  These  equations  are  specialized  for  a  toroidal  shell 
loaded  by  uniform  external  pressure.  The  resulting  equations  are  solved  by  use 
of  series  expansions  in  the  circumferential  and  meridional  directions  for  the 
displacement  components  that  develop  during  buckling.  The  analysis  shows  that 
a  toroidal  shell  can  buckle  in  a  mode  which  is  either  symmetric  or  antimetric 
about  the  equatorial  plane  and  that  the  corresponding  buckling  pressures  are 
always  close  to  each  other.  Axially  symmetric  as  well  as  asymmetric  buckling 
modes  are  considered.  The  numerical  results  show  that  the  asymmetric  modes 
give  lower  buckling  pressures.  Design  curves  which  give  nondimensional  buck¬ 
ling  pressures  for  a  wide  range  of  the  toroidal  shell's  geometric  parameters 
are  presented.  In  addition ,  the  var  iation  of  the  mode  shapes  with  the  geometric 
parameters  is  illustrated.  In  a  comparison  between  the  results  of  the  present 
theory  and  tiv~  few  available  tests  on  toroidal  shells ,  it  is  shown  that  test  and 
theory  agree  to  within  10%.  For  the  limiting  case  of  a  sphere  under  external 
pressure,  the  well-known  classical  result  (p  =  1.21  Eh  /a*)  is  obtained 
numerically  for  both  the  asymmetric  and  axially  symmetric  buckling  modes. 
Finally ,  equations  governing  the  free  vibrations  of  a  prestressed  shell  of  revo¬ 
lution  are  presented  and  specialized  for  a  toroidal  shell  subject  to  external  or 
internal  preseure. 
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